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Executive  Summary 


During  the  course  of  this  project,  our  efforts  have  focused  on  implementing  seismic  event 
characterization  capabilities  at  the  prototype  International  Data  Centre  (PIDC)  and  on  quantifying 
the  performance  of  these  capabilities  with  regard  to  monitoring  the  Comprehensive  Nuclear  Test- 
Ban  Treaty  (CTBT).  We  have  also  been  responsible  for  developing,  operating  and  maintaining  the 
PIDC  World  Wide  Web  site.  Key  accomplishments  of  these  efforts  are  summarized  here  and 
described  in  more  detail  in  the  body  of  this  report. 

First,  we  have  developed  seismic  event-screening  procedures  and  criteria,  based  on  standard  event 
characterization  parameters  computed  at  the  PIDC,  with  the  objective  of  screening  out  events  that 
are  consistent  with  natural  phenomena.  The  screening  criteria  are  based  on  location,  focal  depth, 
and  the  difference  of  body  (mb)  and  surface  (Ms)  wave  magnitudes.  Uncertainties  in  the  estimates 
are  included  in  the  analyses  so  that  events  may  be  screened  out  at  a  specified  confidence  level. 
Location  error  ellipses  are  being  used  on  an  interim  basis  to  categorize  offshore  events  until 
hydroacoustic  screening  capabilities  are  further  developed  and  additional  IMS  hydroacoustic 
stations  are  on-line.  A  standard  set  of  screening  criteria  are  applied  on  a  routine  basis  to  all  events 
above  mb  3.5  in  the  Reviewed  Event  Bulletin  (REB).  In  Section  1  we  describe  the  event 
characterization  parameters,  the  screening  criteria,  and  numerical  indication  of  the  results. 

We  have  also  implemented  and  tested  capabilities  at  the  PIDC  to  compare  events  to  respective 
regional  event  populations  using  high-frequency  (e.g.,  >  4  Hz)  regional  amplitude  ratios,  Pn/Lg 
and  Pn/Sn.  Additional  work  is  needed  to  further  quantify  region-specific  distance  dependence  and 
other  subregional  effects  before  this  analysis  can  be  made  operational.  We  describe  the  procedure 
in  Section  1  and  provide  results  of  applications  to  some  interesting  events  in  Sections  4  and  5. 

Second,  we  have  assessed  the  performance  of  the  screening  criteria  for  all  events  above  mb  3.5  in 
the  REB  for  1996  and  1997.  During  this  two-year  period  there  were  over  30,000  such  events, 
including  two  known  underground  nuclear  tests  at  the  Lop  Nor  test  site  in  China.  We  have  also 
applied  the  screening  criteria  to  45  and  48  historical  underground  nuclear  explosions  at  the  Nevada 
Test  Site  and  Semipalatinsk,  respectively,  using  data  provided  to  us  by  Murphy  (1998).  In  addition, 
we  have  made  comparisons  of  the  results  using  alternative  screening  criteria.  The  results, 
presented  in  Section  2,  indicate  that  a  significant  percentage  of  events  (e.g.,  about  77%  of  the 
events  above  mb  3.5)  can  be  screened  out  by  the  rather  conservative  criteria  considered  and  none 
of  the  95  known  explosions  were  screened  out.  Significant  improvements  in  the  percentage  of 
events  that  can  be  screened  out  are  expected  once  the  screening  analyses  based  on  regional  seismic 
data  are  more  fully  implemented  and  as  the  IMS  seismic  network  is  completed. 
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Third,  we  have  implemented  interactive  capabilities  for  custom  event-screening  that  are  available 
on  the  PIDC  World  Wide  Web  site.  These  on-line  capabilities  allow  remote  users  to  specify  and 
apply  their  own  event  selection  and  screening  criteria,  as  well  as  view  and  download  the  results  of 
the  analyses.  These  capabilities  and  the  data  products  are  described  by  Fisk  et  al.  (1996a).  Carlson 
and  Fisk  (1998)  provide  a  user’s  guide  for  these  capabilities. 

Fourth,  we  have  quantified  the  uncertainties  of  magnitude  (mb  and  Ms)  estimates,  including  the 
dependence  of  the  uncertainties  on  the  number  and  the  coverage  of  stations  used  in  obtaining  the 
magnitude  estimates.  This  work  is  described  by  Bottone  et  al.  (1996). 

Fifth,  we  have  estimated  preliminary  region-specific  distance  corrections  for  regional  seismic 
amplitude  ratios,  Pn/Lg  and  Pn/Sn  in  the  2-4,  4-6  and  6-8  Hz  bands,  for  the  existing  IMS  seismic 
stations.  Results  of  our  initial  work  are  described  by  Bottone  et  al.  (1997).  More  recent  results  of 
progressive  improvements  to  the  distance  corrections,  using  additional  regional  seismic  data,  a 
parametrization  of  the  distance  dependence  which  includes  geometrical  spreading  and  attenuation 
terms,  and  treatment  of  subregional  effects,  are  described  in  Section  3  of  this  report. 

Last,  we  have  performed  in-depth  regional  analysis  of  some  interesting  events  in  the  REB.  These 
events  include  the  two  underground  nuclear  explosions  on  8  June  and  29  July  1996  at  the  Lop  Nor 
test  site  in  China,  an  event  on  10  November  1996  in  Southern  India,  and  an  event  in  the  Kara  Sea 
on  16  August  1997.  The  Lop  Nor  events  are  of  interest  because  they  provide  rare  cases  with  which 
to  test  the  effectiveness  of  the  seismic  event-screening  criteria  with  regard  to  not  screening  out 
simple  underground  explosions.  The  96/07/29  Lop  Nor  explosion  and  the  other  events  are  of 
further  interest  because  the  depth  and  Ms:mb  screening  criteria  cannot  be  applied,  based  on  the 
available  data;  thus  illustrating  the  need  and  utility  of  the  regional  analysis.  In  Section  4  we  present 
results  of  event  characterization  analyses  for  the  Lop  Nor  explosions  and  the  event  in  Southern 
India.  In  Section  5  we  present  results  of  the  analyses  for  the  97/08/16  Kara  Sea  event.  (Note:  This 
type  of  in-depth  event  analysis  is  not  necessarily  considered  within  the  scope  of  the  IDC  but, 
rather,  was  performed  as  a  contractor  to  the  U.S.  Department  of  Defense.) 
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7,603  regional  events.  Superimposed  are  the  scatterplots  as  in  Figure  10  and  the  binned 
means  as  in  Figure  11.  The  table  gives  the  number  of  measurements  for  each  amplitude 
ratio,  the  estimated  standard  deviation  of  the  data  before  and  after  distance  corrections,  and 
the  best-fit  coefficients . 32 
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XIV 


Section  1 

Event-Screening  Approach 


In  this  section  we  describe  the  event  characterization  parameters  that  are  currently  being  computed 
on  a  routine  basis  at  the  PIDC  for  seismic  events  in  the  Reviewed  Event  Bulletin  (REB)  We  also 
describe  the  initial  version  of  the  standard  seismic  event-screening  procedure,  including 
definitions  of  the  event  categories,  screening  criteria,  and  numerical  indication  of  the  results. 

1.1  Event-Characterization  Parameters. 

Here  we  provide  the  definitions  of  seismic  event-characterization  parameters,  in  accordance  with 
Annex  2  of  the  CTBT  Protocol,  that  are  currently  being  computed  and  tested  at  the  PIDC.  (The 
processing  capabilities  to  compute  these  parameters  were  implemented  by  other  contractors.) 

1.1.1  Focal  Depth  and  Epicentral  Location. 

Hypocentral  inversion  is  performed  as  an  iterative  non-linear  least-squares  inversion  of  travel  time, 
azimuth  and/or  velocity  measurements  (see  GSE/CRP/243,  Vol.  2,  Appendix  E).  The  event 
location  and  depth,  confidence  bounds,  and  residuals  are  estimated  using  these  measurements  from 
stations  at  regional  and  teleseismic  distances.  If  the  depth  becomes  negative  for  several  iterations, 
the  depth  estimate  is  fixed  at  zero.  (The  analyst  can  also  constrain  the  depth  estimate  to  the 
surface.)  The  LASPEI91  travel  time  curves  are  used  as  the  global  standard  travel  time  curves  for 
event  location.  Where  available,  regional  travel  time  curves  can  be  incorporated  into  the  analysis. 

Error  estimation  involves  determining  event  error  ellipsoid  and  normalized  confidence  regions, 
utilizing  the  procedure  of  Jordan  and  Sverdrup  (1981),  where  the  error  analysis  can  be  represented 
as  a  combination  of  a  priori  and  a  posteriori  distributions,  through  control  of  two  input 
parameters,  namely,  the  number  of  degrees  of  freedom,  num_dof,  and  an  a  priori  scale  variance 
factor,  est_std_err  (Bratt  and  Bache,  1988).  Given  the  covariance  matrix  of  the  hypocentral 
estimate,  the  error  ellipsoid  and  confidence  regions  are  calculated  from  the  appropriate  marginal 
variances,  based  on  the  F-distribution.  The  default  uses  est_std_err  of  1.0  and  num_ dof  equal  to 
infinity.  In  this  limit,  the  F-distribution  becomes  a  chi-squared  distribution.  The  default  confidence 
level  is  0.90,  which  may  be  modified  to  other  confidence  levels  between  0.0  and  1.0.  The  following 
parameters  are  estimated: 

•  The  two-dimensional  confidence  error  ellipse  of  the  epicenter,  as  determined  from  the 
2x2  covariance  matrix  of  the  epicenter  (marginal  with  respect  to  origin  time  and  depth). 
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•  The  one-dimensional  confidence  interval  of  focal  depth,  as  determined  from  the  scalar 
variance  of  the  focal  depth  (marginal  with  respect  to  origin  time  and  epicenter). 

1.1.2  Body  Wave  Magnitudes  (mb). 

Body  wave  magnitudes,  mb,  are  computed  from  amplitude  and  period  measurements  of  associated 
P  phases  (GSE/CRP/243,  Vol.  2,  Appendix  F).  The  network-average  magnitude  is  based  on  station 
magnitudes  for  each  station  with  the  first  associated  arrival  either  P,  Pn,  or  Pg;  and  it  must  have 
unclipped  amplitude  and  period  data.  Station  magnitudes  are  calculated  with  the  standard  formula: 

mb  =  log Xq(A/T)  +  B(D,  z)  +  S ,  (1) 

where  A  is  the  amplitude,  T  is  the  period,  S  is  a  station  correction  term,  and  the  distance 
corrections,  B(D,  z)  of  Veith  and  Clawson  (1972)  are  used.  The  network-average  magnitude  is  the 
weighted  mean  of  the  station  magnitudes.  Weights  are  inversely  proportional  to  the  sum  of  the 
variances  for  the  station  magnitudes,  station  noise,  and  station  corrections.  Station  corrections  and 
noise  are  stored  in  the  siteaux  table  of  the  PIDC  database.  Error  in  the  network-average  mb 
estimate  is  computed  assuming  a  normal  distribution  of  the  station  magnitude  estimates. 

1.1.3  Surface  Wave  Magnitudes  (Ms). 

The  computation  of  surface  wave  magnitudes,  Ms,  are  carried  out  in  a  manner  similar  to  that  used 
for  the  mb  calculations.  The  error  is  also  estimated  in  a  similar  way.  The  formula  for  Ms  station 
magnitudes  is  the  “Prague  formula”  of  Vanek  et  al.  (1962): 

Ms  =  log \o(A/T)  +  1.66  log10(D)  +  0.3  .  (2) 

where  A  is  the  instrument-corrected  zero-to-peak  amplitude  (in  nanometers)  with  period,  T, 
between  18  and  22  seconds,  and  D  is  the  distance.  The  arrival  is  also  required  to  be  within  a 
velocity  window  of  ±0.3  km/sec  of  the  predicted  arrival  time  based  on  a  global  dispersion  model. 

1.1.4  Teleseismic  P-Wave  Complexity. 

The  measure  of  complexity  implemented  at  the  PIDC  follows  the  energy-ratio  technique  of 
Douglas  (1980).  In  this  method  the  coherent  beam  of  the  P  signal  is  filtered  in  the  band  from  0.8 
to  2.0  Hz,  squared  and  smoothed  with  an  exponential  window  (with  a  time  constant  of  1.5  sec). 
Complexity  is  then  computed  as  the  ratio  of  the  area  under  the  above  curve  in  the  signal_len2 
seconds  after  the  P  signal  to  the  area  in  the  signaljenl  seconds  of  the  P  signal  ( S ),  both  corrected 
for  background  noise  {N).  Under  the  current  default,  signaljenl  =  5  and  signal Jen2  =  30,  and 
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(3) 


Complexity  = 


5(5  -  35  seconds)  -  N 
5(0  -  5  seconds)  -  N 


with  the  exception  of  the  following  two  cases: 


•  If  5(0-5  sec)  -  N  <  0 ,  complexity  is  undefined  and  so  is  set  to  -1; 

•  If  5(5-55  sec)  -  N  <  0 ,  all  energy  is  in  the  first  5  sec,  so  complexity  is  set  to  0.0. 


Complexity  is  measured  only  for  stations  that  are  associated  to  an  event  and  are  at  teleseismic 
distances  (30-90  degrees).  If  a  PcP  phase  is  observed  within  35  sec  of  P,  no  measurement  is  made. 
The  complexity  measurements  are  stored  in  the  complexity  table  of  the  PIDC  database. 


1.1.5  Regional  Phase  Time-Domain  Amplitudes. 

The  Detection  and  Feature  Extraction  (DFX)  code  at  the  PIDC  can  compute  time-domain 
amplitude  measurements  on  a  variety  of  beams  for  specified  time  windows.  Amplitude 
measurement  types  include  absolute  maximum  amplitude,  maximum  peak-to-trough,  maximum 
peak-to-peak,  RMS  amplitude,  etc.  Beam  types  include  coherent,  steered,  incoherent,  and  RMS 
measured  on  vertical,  radial,  or  transverse  components.  Time  intervals  for  the  measurements  can 
be  based  on  predicted  travel  times  of  seismic  phases,  or  on  fixed  group  velocity  windows,  or  they 
may  be  computed  from  observed  arrival  times  in  the  database. 

Currently,  the  following  time-domain  regional  phase  amplitude  measurements  are  computed  for 
each  station  that  is  associated  to  and  within  20  degrees  of  an  event: 

a)  Absolute  maximum  amplitude  on  2-4,  4-6,  6-8,  8-10  Hz  RMS  beams  for  predicted  time  and/or 
velocity  windows  around  Pn,  Pg,  Sn,  Lg  phases: 

Pn:  8s  before  the  theoretical  arrival  time  of  Pn  to  a  group  velocity  of  6.4  km/s 

Pg:  6.3  to  5.8  km/s  group  velocity 

Sn:  5s  before  the  theoretical  arrival  time  of  Sn  and  a  20s  duration 

Lg:  3.7  to  3.0  km/s  group  velocity 

b)  Absolute  maximum  amplitude  on  2-4, 4-6,  6-8,  8-10  Hz  RMS  beams  for  predicted  time  and/or 
velocity  windows  around  pre-Pn,  pre-Pg,  pre-Sn,  and  pre-Lg  noise: 

Pre-Pn:  13s  before  the  theoretical  arrival  time  of  Pn  and  a  5s  duration 

Pre-Pg:  6.4  to  6.3  km/s  group  velocity 

Pre-Sn:  10s  before  the  theoretical  arrival  time  of  Sn  and  a  5s  duration 

Pre-Lg:  3.8  to  3.7  km/s  group  velocity 
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The  instrument  response  correction  at  the  center  frequency  is  applied  to  all  of  the  amplitude 
measurements.  Also,  amplitudes  in  the  8-10  Hz  band  are  only  calculated  for  stations  with  a 
Nyquist  frequency  >  10  Hz.  These  signal  and  noise  amplitude  measurements  are  stored  in  the 
originamp  table  of  the  PIDC  database. 

1.1.6  Short-Period/Long-Period  Energy  Ratio  (SP/LP). 

The  SP/LP  event-characterization  parameter  (e.g.,  Saikia  et  al.,  1996)  is  the  ratio  of  the  integrated 
energy  between  the  short-period  (SP)  vertical  component  P-wave  train  and  the  sum  of  the  three 
long-period  (LP)  surface  wave  trains.  The  calculation  of  the  energy  ratio  is  restricted  to  3- 
component  broadband  stations  within  15  degrees  of  the  event.  The  P-wave  train  is  the  filtered  (1-8 
Hz)  time  series  between  a  group  velocity  of  8.2  km/s  and  4.5  km/s,  while  the  LP  surface  wave 
trains  are  the  filtered  (0.033  -  0.167  Hz)  Z,  N  and  E  channel  time  series  between  a  group  velocity 
of  4.5  km/s  and  2.5  km/s.  SP/LP  energy  ratio  values  can  be  found  in  database  table  splp. 

1.1.7  Spectral  and  Cepstral  Measurements. 

Spectral  and  cepstral  analyses  are  performed  on  each  associated  regional  phase  within  20  degrees 
(Baumgardt  and  Zeigler,  1988).  The  following  measurements  are  made  and  stored  in  the  database: 

(a)  The  variance  of  the  detrended  log  spectrum  of  each  associated  regional  phase  is  calculated  for 
frequency  bands  with  width  4.0  Hz  and  with  SNR  3.0.  The  variance  and  the  coefficients  of  the 
quadratic  trend,  removed  from  the  log  spectrum,  are  stored  in  table  spvar. 

(b)  Amplitudes  and  quefrencies  of  all  cepstral  peaks  with  amplitudes  above  a  default  threshold  of 
0.032,  that  are  not  present  in  the  preceding  noise  window,  are  calculated.  A  5.0  sec  window  starting 
0.3  sec  before  the  time  of  the  associated  phase  is  used  for  the  signal,  and  the  preceding  5.0  sec 
window  is  used  for  the  noise.  Cepstral  peaks  are  only  declared  if  they  are  consistent  (i.e.,  the  same 
quefrency)  among  two  or  more  phases  at  the  same  station.  The  maximum  amplitude  of  the 
individual  phase  peaks  for  that  station  and  the  quefrency  are  stored  in  database  table  ceppks. 

1.2  Screening  Methodologies  for  Seismic  Events. 

Seismic  event-screening  criteria  that  are  currently  applied  on  a  routine  basis  to  the  events  in  the 
REB  are  based  on  depth,  Ms:mb,  and  location  (to  determine  offshore  events).  Here  we  describe  the 
screening  methodologies  based  on  these  measurements.  We  also  briefly  describe  the  approach  to 
compare  events  to  regional  event  populations,  based  on  Pn/Lg  and  Pn/Sn. 
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1.2.1  Screening  Methodology  based  on  Depth. 

Given  the  depth  estimate,  D ,  and  its  variance,  aD ,  ( depth  and  szz  in  the  database  schema)  from  the 
hypocentral  inversion  described  above,  the  distribution  of  the  square  of  the  depth  estimate 
(normalized  to  zero  mean  and  unit  variance)  is  obtained  as  the  marginal  distribution  of  the 
hypocenter  estimate  with  respect  to  origin  time  and  epicenter.  This  is  a  chi-squared  distribution 
with  one  degree  of  freedom,  which  is  equivalent  to  a  normal  distribution  for  the  depth  estimate. 
That  is,  (D  -  D)/cd  has  a  normal  distribution  with  zero  mean  and  unit  variance,  where  D  is  the 
true  (unknown)  depth.  Under  this  assumption,  a  one-sided  100(1  —  oc)%  confidence  interval  for 
the  depth,  D,  is  [D  -  za <JD,  00 ) ,  where  za  is  the  a-percentile  of  the  normal  distribution  with  zero 
mean  and  unit  variance.  That  is,  za  is  defined  such  that  P[z  >  zj  =  l-<P(z0)  =  a ,  where 

Z  00 

<b(z)  =  -L  f  dxt*/2  =  1  -  -L \dxt~x2/2.  (4) 

V2iJ 

-00  z 

An  event  is  screened  out  at  the  100(  1  -  a)%  confidence  level  if  the  shallowest  part  of  the  depth 
confidence  interval  is  deeper  than  the  depth  threshold,  i.e.,  if 

b-zaoD>DQ.  (5) 

The  default  depth  threshold  is  set  at  10  km.  For  a  confidence  level  of  97.5%,  za  =  1.96 . 


1.2.2  Screening  Methodology  based  on  Ms:mb. 


For  an  event  for  which  there  are  both  mb  and  Ms  magnitude  measurements,  let  Nb  be  the  number 
of  stations  ( nsta  in  the  database  schema)  with  mb  measurements.  The  network  average  of  mb,  mb , 
is  given  by  the  sample  mean  of  the  Nb  individual  station  magnitude  measurements,  mb  , 

1 

mb=TFyLmb,i’  (6) 

b  i=  1 

The  network  average  of  Ms,  Ms ,  is  computed  in  a  slightly  more  complicated  manner.  For  a  given 
station,  j ,  let  n-  be  the  number  of  elements  at  a  station  which  measured  Ms  and  let  Ns  be  the  total 
number  of  stations  measuring  Ms.  The  station  average,  M s  • ,  is  given  by  the  sample  mean  of  the 
n  j  values  of  Ms  at  station  j ,  and  the  network  average  of  Ms  is  given  by 


The  network-average  estimate  of  Amb-Ms  is  then  taken  to  be  Amb  -  Ms ,  where  A  is  a  constant 
slope  term  to  account  for  possible  magnitude  dependence  of  the  Ms:mb  relation. 
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The  standard  deviations  for  the  individual  mb  and  Ms  measurements  are  currently  taken  to  be 

ob  =  <3 s  =  0.3  ,  independent  of  the  value  of  mb  or  Ms.  This  is  an  upper  bound  of  the  standard 

deviations  estimated  by  Bottone  et  al.  (1996).  The  variance  of  the  mean  of  Nb  independent 

2 

measurements  of  mb  is  given  by  cb/Nb,  with  a  similar  expression  for  Ms.  The  variance  of 

—  2 

Amb  -  Ms ,  denoted  by  c~M ,  assuming  that  the  mb  and  Ms  measurements  are  uncorrelated,  is 
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varian ce(A?nb  -  M s) 


(8) 


_  —  2 
Under  the  assumption  that  Amb-Ms  has  a  normal  distribution  with  variance  gm,  a  one-sided 

100(1  -  a)%  confidence  interval  for  Amb-Ms  is  (-°°,  Ainb-Ms  +  za<3M].  An  event  is  screened 

out  if  this  confidence  interval  is  entirely  less  than  a  threshold,  M0,  i.e.,  if  it  is  entirely  outside  of 

the  explosion  population.  Thus,  an  event  is  screened  out  if 


Amb-Ms  +  zaoM<M0.  (9) 

Marshall  and  Basham  (1972)  have  suggested  a  threshold  of  mb-Ms  =  1.2,  based  on  historical 
explosion  and  earthquake  data.  After  careful  calibration  of  PIDC  magnitudes,  relative  to  historical 
data,  Murphy  (1997)  recommends  that  a  threshold  of  1.25mb-Ms  =  2.20  be  used  at  the  PIDC. 
Testing  of  both  of  these  criteria,  using  a  97.5%  confidence  level,  is  described  in  Section  2. 


1.2.3  Numerical  Indication  of  the  Screening  Results. 

The  CTBT  Protocol  states  that,  “The  standard  event  bulletin  shall  indicate  numerically  for  each 
event  the  degree  to  which  that  event  meets  or  does  not  meet  the  event-screening  criteria.”  In  light 
of  this  requirement,  probabilities  are  computed  to  indicate  the  degree  to  which  the  depth  and 
Ms:mb  estimates  are  inconsistent  with  possible  values  from  a  shallow  explosion  source, 
accounting  for  the  uncertainties  in  the  estimates.  An  alternative  representation  of  the  probabilities 
in  terms  of  a  screening  “score”  is  also  presented. 

The  screening  criterion  for  depth  in  Equation  (5),  in  terms  of  the  depth  confidence  interval,  is 
equivalent  to  performing  a  test  of  hypothesis.  That  is,  under  the  null  hypothesis  that  the  true  depth, 
D,  is  less  than  or  equal  to  a  threshold  depth,  D0,  the  probability  of  obtaining  an  estimate,  D ,  deeper 
than  some  fixed  depth,  D,  given  the  standard  deviation,  oD,  is  P[D>  D\D0,gd]  =  l-d>(x), 
where  x  =  (D  -  Dq)/gd  and  <!>  ( • )  is  defined  in  Equation  (4).  Thus,  if  the  measured  value  of  the 
depth  is  Dx ,  the  hypothesis  that  the  true  depth  of  the  event  was  less  than  or  equal  to  D0  may  be 
rejected  at  the  a  significance  level  if  P[D  >  dAd0,  Od ]  <  a  or,  equivalently,  if 
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Pdepth^nb<Dl\D0,oD]  =  «&(*)>! -ot. 


(10) 


Note  that  P[D  >  Dx  Jz>0,  <5D\  is  referred  to  as  the  “p- value”  in  the  statistics  literature  and  Pdepth 
is  its  complement.  An  equivalent  “score”  for  the  depth  screening  criterion  may  also  be  defined  as 

SC0RE depth  =  (A  -  t>oVZa?D  ~  1  =  ^cc  "  1  •  d  D 

Note  that  a  positive  score  corresponds  to  Pdepth  >  1  -  a  and  a  one-sided  100(  1  -  a)%  confidence 
interval  for  depth  entirely  deeper  than  D0.  For  a  one-sided  97.5%  confidence  interval,  za  =  1.96 . 
Thus,  the  score  increments  by  one  for  an  increment  of  roughly  two-sigma  from  the  decision  line. 

In  this  analysis  we  assume  that  the  errors  are  distributed  according  to  a  Gaussian  (normal) 
distribution  for  positive  depths,  with  zero  probability  for  values  above  the  ground.  (Recall  from 
Section  1.1.1  that  negative  depth  estimates  are  constrained  to  the  surface.)  We  also  assume  that  the 
variance  is  independent  of  depth  and  that  the  value  szz  used  for  oD  sufficiently  represents  the  true 
uncertainty.  In  the  future  we  hope  to  verify  or,  if  necessary,  improve  on  these  assumptions. 

Similarly,  the  probability  of  obtaining  an  estimate  of  Amb-Ms,  M ,  less  than  some  fixed  value,  M , 
given  that  the  true  value  is  at  the  threshold  value,  M0,  and  that  the  standard  deviation  is  oM ,  is 
given  by  P[M  <  M\M0,  oM]  =  1  -  3>(y) ,  where  y  =  (M0  -  M)/ oM .  Thus,  if  a  measurement  of 
Amb-Ms  is  M\ ,  the  null  hypothesis  that  the  true  value  of  Amb-Ms  is  greater  than  or  equal  to  M0 
may  be  rejected  at  the  a  significance  level  if  P[M  <  M\\M0,  aM]  <  a  or,  equivalently,  if 

Rmbms  ~  P[M  >  Mi  |M0,  oM]  =  O(y)  >  1  -  a .  (12) 

An  equivalent  “score”  for  the  Ms:mb  screening  criterion  may  also  be  defined  as 

SCOREmbms  =  (Mq  -  Mj  )/zacM  -  1  =  y/za  -  1  •  (13) 

Assuming  the  depth  and  Ms:mb  estimates  are  approximately  independent  (which  is  valid  under  the 
null  hypothesis  that  the  event  is  shallow),  the  combined  probability  that  a  shallow  (D  <  D0) 
explosion  would  have  given  a  depth  measurement  greater  than  or  equal  to  Dx  and  a  measurement 
of  Amb-Ms  less  than  Mx  is  given  by  P  =  P[D  >  Dx  |d0,  aD]  ■  P[M  <  M] \M0,  <jm]  .  (It  is 
straightforward  to  modify  this  expression  to  treat  the  correlation  between  the  depth  and  Ms:mb,  if 
the  correlation  is  quantified.)  This  relation  can  be  used  to  test  the  hypothesis  that  the  event  could 
be  a  shallow  explosion,  based  on  depth  and  Ms:mb.  If  the  hypothesis  is  rejected,  the  event  would 
be  screened  out.  The  complement  of  this  relation  is  given  by 

R combined  =  ^  ~  ^  ~  R depth>^  ~  R mbtns^  ■  (14) 
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An  equivalent  “score”  for  the  combined  depth  and  Ms:mb  screening  criteria  may  also  be  defined  as 


SCORE combined  =  z/za  -  1 .  (15) 

where  z  is  determined  from  the  expression  Pcombined  =  O(z) . 

The  expressions  for  P depth ,  P mbms ,  and  Pcombined  are  use^  t0  numerically  indicate  the  screening 
results.  Each  can  have  possible  values  between  0  and  1,  with  a  value  of  0  corresponding  to  the 
measurements  entirely  within  the  possible  explosion  range  and  a  value  of  1  corresponding  to  the 
measurements  entirely  outside  of  the  possible  explosion  range,  accounting  for  the  uncertainties. 
These  expressions  can  be  generalized  to  an  arbitrary  set  of  event-screening  criteria. 

1.2.4  Screening  Methodology  based  on  Location. 

Location  confidence  ellipses  are  being  used  on  an  interim  basis  to  categorize  offshore  events  until 
additional  IMS  hydroacoustic  stations  are  on-line  and  hydroacoustic  event-screening  capabilities 
are  further  developed  at  the  PIDC.  An  onshore/offshore  grid  (including  islands)  with  one-degree 
resolution  is  currently  used.  The  90%  location  confidence  ellipse,  obtained  from  the  hypocentral 
inversion,  must  not  touch  or  contain  any  “onshore”  cells  for  an  event  to  be  categorized  as  offshore. 
The  percentage  of  offshore  grid  cells  enclosed  or  touching  the  ellipse,  pct_offshore,  is  computed. 
The  grid,  confidence  ellipses,  and  results  of  applications  are  illustrated  in  Section  2.2. 

1.2.5  Screening  Methodology  for  Regional  Seismic  Events. 

The  screening  methodology  for  regional  seismic  events  is  based  on  the  likelihood-ratio  outlier  test, 
as  described  in  detail  by  Fisk  et  al.  (1996b).  The  analysis  utilizes  high-frequency  Pn/Lg  and  Pn/Sn 
measurements,  described  above,  to  compare  events  to  the  relevant  regional  event  population.  The 
p-value  of  the  hypothesis  test  is  used  to  numerically  indicate  the  screening  result.  The 
methodology  has  been  implemented  and  tested  on  a  developmental  basis  at  the  PIDC  (e.g.,  Fisk  et 
al.,  1996a).  Further  work  is  needed  to  establish  regional  distance  corrections  and  grids  of  regional 
phase  blockage  before  this  analysis  is  ready  to  be  integrated  into  operations  at  the  PIDC. 

1.3  Initial  Version  of  the  Seismic  Event-Screening  Procedure. 

We  now  describe  the  initial  version  of  the  standard  event-screening  procedure  that  was  applied  to 
all  REB  events  on  a  routine  daily  basis  during  1997,  including  definitions  of  the  event  categories, 
the  explicit  screening  criteria,  the  input  data  and  parameters,  and  the  output. 
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1.3.1  Category  Definitions  and  Event-Screening  Criteria. 

Table  1  defines  the  categories  and  experimental  event-screening  criteria,  based  on  location,  depth 
and  Ms:mb,  that  were  applied  on  a  routine  basis  to  all  events  in  the  REB  for  1997.  Events  below 
mb  3.5  are  not  considered  for  the  screening  analyses  because  the  screening  criteria  that  are 
currently  used  are  typically  applicable  to  events  that  are  large  enough  so  that  depth  and  mb-Ms 
measurements  can  be  obtained,  and  because  the  existing  screening  criteria  do  not  include 
capabilities  to  screen  out  chemical  mining  blasts,  which  are  typically  smaller  than  mb  3.5. 

Events  with  90%  location  confidence  ellipses  that  are  not  entirely  offshore,  and  for  which  mb-Ms 
and  unconstrained  depth  estimates  are  not  available,  fall  in  the  insufficient  data  category,  a  subset 
of  those  events  that  are  not  screened  out.  Events  with  location  confidence  ellipses  that  are  entirely 
offshore  and  with  constrained  depths  and  no  mb-Ms  are  placed  in  a  separate  category,  “Offshore 
Only”.  These  could  be  considered  as  offshore  events  with  insufficient  data. 


Table  1.  Seismic  event-screening  categories  and  criteria. 


Category 

Criteria 

Not  Considered 

•  mb  <  3.5 

Insufficient  Data 

•  90%  location  error  ellipse  not  entirely  offshore  and 

•  depth  constrained  to  the  surface  and  no  mb-Ms 

Offshore  Only 

•  90%  location  error  ellipse  entirely  offshore  and 

•  depth  constrained  to  the  surface  and  no  mb-Ms 

Screened  Out 

•  97.5%  confidence  interval  for  depth  >  10  km  or 

•  97.5%  confidence  interval  for  mb-Ms  <  1.2 

Not  Screened  Out 

•  90%  location  error  ellipse  is  not  entirely  offshore  and 

•  97.5%  confidence  interval  for  depth  not  >  10  km  and 

•  97.5%  confidence  interval  for  mb-Ms  not  <1.2 

Note  that  the  one-sided  97.5%  confidence  intervals  are  equivalent  to  the  two-sided  95%  confidence 
intervals  that  have  been  utilized  and  reported  on  previously.  Also,  as  of  1  January  1998,  the  Ms:mb 
screening  criterion  is  1.25mb-Ms  <  2.20,  as  recommended  by  Murphy  (1997).  Other  changes  to 
these  provisional  event-screening  criteria  may  also  be  made  in  the  future,  based  on 
recommendations  by  Working  Group  B  of  the  CTBT  Preparatory  Commission. 
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1.3.2  Input  Data  and  Parameters. 


Table  2  provides  a  list  of  the  data  used  in  the  existing  version  of  the  event-screening  code,  as  well 
as  the  database  tables  in  which  they  are  stored  and  their  descriptions.  These  data  are  obtained  from 
the  IDCREB  database  account,  after  the  seismic  analysts  have  completed  their  review  of  the  events. 

Table  3  lists  the  input  parameters  to  this  version  of  the  code,  their  default  values  (used  during 
1997)  and  descriptions.  The  parameter,  Idpflg,  can  be  set  to  true  to  require  that  one  or  more  depth 
phases  be  detected  for  the  depth  screening  analysis  to  be  applied.  The  parameter,  nstaO,  is  the 
minimum  number  of  stations  required  to  have  Ms  measurements  for  the  Ms:mb  screening  criterion 
to  be  applied.  The  user  may  also  specify  a  custom  rectangular  region  and  time  interval  of  interest. 


Table  2.  Input  data  to  the  event-screening  procedure. 


Attribute  Name 

Database  Table 

Attribute  Description 

lat 

origin 

estimated  latitude 

Ion 

origin 

estimated  longitude 

depth 

origin 

estimated  depth 

time 

origin 

epoch  origin  time 

orid 

origin,  origerr 

origin  identification  number 

ndp 

origin 

number  of  depth  phases  detected 

mb 

origin 

body  wave  magnitude 

mbid 

origin,  netmag 

mb  magnitude  identifier  (magid) 

ms 

origin 

surface  wave  magnitude 

msid 

origin,  netmag 

ms  magnitude  identifier  (magid) 

smajax 

origerr 

semi-major  axis  of  location  error 

sminax 

origerr 

semi-minor  axis  of  location  error 

strike 

origerr 

strike  of  semi-major  axis 

szz 

origerr 

variance  of  depth  estimate 

nsta 

netmag 

number  of  stations  used  in  magnitude  estimates 
(corresponding  to  magid  =  mbid  and  magid  =  msid) 
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Table  3.  Input  parameters  to  the  event-screening  procedure. 


Parameter 

Default  Value 

Description 

confloc 

0.90 

confidence  level  for  the  location  error  ellipse 

confdep 

0.975 

confidence  level  for  the  depth  criterion 

confmbms 

0.975 

confidence  level  for  the  Ms:mb  criterion 

depthO 

10.0 

depth  threshold  in  km 

ambmsO 

1.2 

threshold  for  (A  mb  -  Ms) 

amb_cal 

1.0 

slope  (A)  of  the  (A  mb  -  Ms)  relation 

mbO 

3.5 

magnitude  (mb)  cut-off  for  considered  events 

Idpflg 

F 

logical  flag  to  require  depth  phases 

nstaO 

1 

Minimum  number  of  stations  required  for  Ms 

1.3.3  Output  Data. 

The  following  computed  quantities  for  each  event  are  the  output  of  the  event-screening  code: 
P depth  P mbms’  P combined’  as  defined  in  Section  1.2.3;  pctjojf shore,  as  defined  in  Section  1.2.4;  and 
the  Category,  as  defined  in  Table  1.  The  numbers  of  events  in  each  category  are  also  compiled  on 
a  daily  basis  and  provided  as  output  to  the  Executive  Summary  (see  Fisk  et  al.,  1996a). 


Section  2 

Results  of  Experimental  Testing  at  the  Prototype  IDC 


2.1  Introduction. 

Here  we  describe  results  of  experimental  testing  of  the  event-screening  capabilities  at  the  PIDC. 
First,  we  apply  a  “standard”  set  of  screening  criteria  to  all  events  in  the  REB  during  1997,  and 
summarize  the  results.  Second,  we  apply  alternative  screening  criteria  for  depth  and  Ms:mb  to  the 
events  in  the  REB  during  1996  and  1997,  and  compare  the  results  to  those  obtained  using  the 
standard  criteria.  During  this  two-year  period  there  were  over  30,000  events  in  the  REB  above  mb 
3.5,  including  two  known  underground  explosions  at  the  Lop  Nor  test  site  in  China.  We  also  apply 
the  screening  criteria  to  45  and  48  historical  nuclear  explosions  at  the  Nevada  Test  Site  and 
Semipalatinsk,  respectively,  using  data  provided  to  us  by  Murphy  (1998).  Third,  we  compute  the 
screening  probabilities  and  scores  (as  defined  in  Section  1.2.3)  for  these  events  and  describe  the 
results.  Last,  we  provide  some  conclusions  and  recommendations. 

2.2  Event-Screening  Performance  Summaries. 

Figure  1  shows  a  map  of  the  existing  IMS  primary  stations  (squares)  and  the  screening  results  for 
20094  REB  events  during  1997.  Event  markers  are  color-coded  by  screening  category  as  defined 
in  the  legend.  Most  of  the  events  below  mb  3.5,  that  were  detected  by  the  existing  IMS  network, 
occurred  in  Scandinavia,  the  U.S.,  western  South  America,  or  Australia,  where  the  network  density 
is  relatively  high.  The  majority  of  events  that  were  not  screened  out  or  had  insufficient  data 
occurred  along  the  Pacific  Rim,  in  Southern  Europe  or  Central  Asia,  where  a  high  percentage  of 
the  events  occur  and,  in  some  cases,  where  the  network  density  is  relatively  sparse  presently. 

Figure  2  summarizes  the  numbers  (upper  plot)  and  percentages  (lower  plot)  of  events  in  each 
screening  category  per  month  for  all  of  1997,  based  on  the  screening  criteria  in  Table  1  (above). 
There  were  typically  between  1400  and  2100  events  in  the  REB  per  month  (black  bars),  based  on 
the  existing  IMS  seismic  network.  About  76%  to  83%  of  these  events  were  mb  >  3.5  (cyan  bars). 

To  illustrate  how  the  results  depend  on  event  magnitude.  Figure  3  shows  the  percentages  of  events 
per  month  by  screening  category  relative  to  the  number  of  events  above  mb  3.5  (upper  plot)  and 
above  mb  4.5  (lower  plot).  The  upper  plot  in  Figure  3  indicates  that,  as  of  1  May  1997,  between 
74%  to  77%  of  the  events  above  mb  3.5  were  screened  out,  roughly  18%  to  23%  have  insufficient 
data  with  which  to  perform  the  screening  analyses,  and  approximately  4%  were  not  screened  out 
based  on  the  experimental  screening  criteria.  The  lower  plot  in  Figure  3  indicates  that,  as  of  1  May 
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1997,  between  88%  to  98%  of  the  events  above  mb  4.5  were  screened  out,  roughly  3%  to  10%  have 
insufficient  data  with  which  to  perform  the  screening  analyses,  and  approximately  3%  or  less  were 
not  screened  out  based  on  the  current  screening  criteria. 

Note  that  three  events  in  the  “insufficient  data”  category  during  1996  and  1997  were  a  96/07/29 
Lop  Nor  explosion,  an  event  in  Southern  India  on  96/1 1/10,  and  an  event  in  the  Kara  Sea  on  97/08/ 
16.  None  of  these  events  were  screened  out;  however,  these  cases  illustrate  the  need  for  screening 
analyses  based  on  regional  data,  as  described  further  in  Sections  4  and  5. 

The  upper  plot  in  Figure  4  shows  more  detailed  dependence  of  the  screening  results  on  mb  for  the 
REB  events  during  JUL-DEC  1997.  The  solid  black  line  depicts  the  total  number  of  REB  events 
during  this  period.  The  cyan  curve  indicates  the  cumulative  number  of  events  above  the  given  mb 
level.  The  solid  blue,  red  and  green  curves  depict  the  cumulative  numbers  of  events  versus  mb  with 
insufficient  data,  that  were  not  screened  out,  and  that  were  screened  out,  respectively.  The  dashed 
and  dotted  curves  indicate  the  cumulative  numbers  of  events  that  were  screened  out  by  the 
individual  depth  and  mb-Ms  criteria,  respectively.  The  intermittent  dotted-dashed  curve  indicates 
the  cumulative  numbers  of  events  that  were  screened  out  by  the  depth  and/or  mb-Ms  criteria. 

In  the  lower  plot  in  Figure  4,  the  cumulative  numbers  have  been  converted  to  percentages,  relative 
to  the  total  number  of  events  above  a  given  magnitude.  This  plot  indicates  that  about  77%  of  the 
events  above  mb  3.5  and  about  91%  of  the  events  above  mb  4.5  were  screened  out.  It  also  depicts 
decreasing  percentages  of  events  that  are  not  screened  out  or  have  insufficient  data  as  functions  of 
increasing  mb.  The  dashed  and  dotted  curves  indicate  that  about  31%  and  24%  of  the  events  above 
mb  3.5  can  be  screened  out  by  depth  and  mb-Ms,  individually,  while  about  46%  of  those  events 
can  be  screened  out  by  depth  and/or  mb-Ms,  with  increasing  percentages  with  increasing  mb.  The 
gap  between  the  total  percentage  of  events  screened  out  and  the  percentage  of  events  screened  out 
by  depth  and/or  mb-Ms  corresponds  to  events  that  were  offshore  at  the  90%  confidence  level. 

Table  4  lists  the  numbers  and  percentages  of  events  screened  out  by  all  combinations  of  individual 
and  combined  depth  and  mb-Ms  screening  criteria  for  the  REB  events  during  1997.  Of  the  20094 
events,  15998  (80%)  were  mb  3.5,  of  which  8017  (50%)  had  mb-Ms  and/or  unconstrained  depth 
estimates  for  which  at  least  one  of  the  depth  and  mb-Ms  screening  criteria  could  be  applied. 
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Table  4.  Numbers  and  percentages  of  1997  REB  events  (above  mb  3.5)  screened  out  by  all 
combinations  of  individual  and  combined  depth  and  mb-Ms  screening  criteria. 


Screening  Criteria 

Number 

Available 

Number 

Screened 

Percentage  of 
Available 

Percentage  of 
mb  3.5 

depth 

5641 

4734 

84% 

30% 

mb-Ms 

4079 

3347 

82% 

21% 

depth  AND  mb-Ms 

1703 

1221 

72% 

7% 

depth  OR  mb-Ms 

8017 

6860 

86% 

43% 

The  results  in  Table  4  indicate  that  using  the  criteria  for  depth  or  for  mb-Ms,  43%  of  all  events 
above  mb  3.5  can  be  screened  out,  while  only  7%  can  be  screened  out  if  events  must  satisfy  both 
the  depth  and  mb-Ms  criteria.  This  result  is  to  be  expected  since  deep  seismic  events  typically  do 
not  have  large,  if  any,  detectable  long-period  surface  waves.  Thus,  it  would  be  very  ineffective,  in 
terms  of  screening  out  events,  to  require  that  both  the  depth  and  mb-Ms  criteria  be  satisfied. 

The  offshore-location  screening  algorithm  (cf.  Section  1.2.4)  has  also  been  tested.  Figure  5 
illustrates  the  grid  and  results  for  events  near  Australia.  The  markers  and  90%  error  ellipses  in 
green  or  red  correspond  to  events  that  are  categorized  as  offshore  or  not  offshore,  respectively. 
Many  of  the  events  with  large  error  ellipses  were  not  screened  out  because  they  touch  or  contain 
“onshore”  grid  cells. 

For  8152  REB  events  (mb  >  3.5 )  that  occurred  between  July  and  December  1997,  4743  events 
(58%)  were  categorized  as  offshore,  and  2238  of  these  events  (27%)  had  constrained  depths  and 
no  Ms.  These  2238  events  are  placed  in  the  “Offshore  Only”  category,  i.e.,  offshore  events  with 
insufficient  data  to  apply  the  seismic  depth  and  Ms: mb  screening  criteria.  In  the  future,  we  plan  to 
develop  hydroacoustic  event-screening  criteria  for  application  to  such  events. 

2.3  Comparison  of  Alternative  Screening  Criteria. 

In  addition  to  the  criteria  in  Table  1  (above),  other  screening  criteria  for  depth  and  Ms:mb  were 
tested  using  the  REB  events  for  1997  with  mb  >  3.5 .  First,  alternative  depth  criteria  were  tested 
using  5641  events  with  unconstrained  depth  estimates.  Requiring  that  the  99%  confidence  interval 
for  depth  be  deeper  than  a  20  km  threshold,  3517  events  are  screened  out,  while  4300  events  are 
screened  out  using  a  99%  confidence  interval  and  a  10  km  threshold.  For  comparison,  4734  events 


14 


are  screened  out  using  a  95%  depth  confidence  interval  and  a  10  km  threshold.  This  illustrates  the 
dependence  of  the  screening  results  on  the  thresholds  and  confidence  levels  that  are  used. 

Figure  6  shows  a  comparison  of  two  different  Ms:mb  criteria  for  all  REB  events  with  mb  >  3.5 
and  Ms  for  1996  (2105  events;  upper  plots)  and  1997  (4079  events;  lower  plots).  The  increase  in 
the  number  of  such  events  for  1997  is  due  to  an  improved  procedure  for  detecting  long-period 
surface  waves  that  was  implemented  early  in  1997  (Stevens,  1997).  The  plots  on  the  left 
correspond  to  a  screening  criterion  of  requiring  that  the  95%  confidence  interval  for  mb-Ms  is 
entirely  less  than  1.2,  while  the  plots  on  the  right  correspond  to  a  95%  confidence  interval  for 
1.25mb-Ms  entirely  less  than  2.20.  The  latter  Ms:mb  criterion  is  based  on  a  careful  calibration 
study  by  Murphy  (1997),  treating  systematic  differences  between  PIDC  and  NEIS  magnitudes  for 
historical  explosions  and  earthquakes,  as  well  as  the  dependence  of  the  Ms:mb  screening  criterion 
on  mb. 

Under  the  two  criteria,  1743  (82.8%)  versus  1629  (77.4%)  events  for  1996,  and  3347  (82.1%) 
versus  3242  (79.5%)  events  for  1997,  were  screened  out,  respectively.  Included  on  both  upper  plots 
of  Figure  6  are  the  Ms  and  mb  values  for  the  96/06/08  Lop  Nor  explosion,  which  is  not  screened 
out  by  either  criterion.  In  fact,  in  both  cases  the  99%  confidence  intervals  are  completely  above  the 
respective  thresholds.  Included  on  both  lower  plots  of  Figure  6  are  the  Ms  and  mb  values  for  45 
and  48  historical  nuclear  explosions  at  NTS  and  Semipalatinsk,  respectively.  The  magnitude  data 
have  been  carefully  calibrated  by  Murphy  (1997)  to  be  consistent  with  PIDC  magnitude  estimates. 
Note  that  none  of  the  explosions  were  screened  out  by  the  latter  Ms:mb  criterion,  while  5  of  93 
explosions  were  screened  out  by  the  former  criterion,  as  depicted  by  the  magenta  squares.  Thus, 
although  slightly  fewer  events  were  screened  out  by  the  latter  Ms:mb  criterion,  the  calibration 
work  by  Murphy  (1997)  and  the  results  presented  here  indicate  that  it  is  more  generally  valid,  in 
terms  of  not  screening  out  explosions,  over  a  wider  range  of  PIDC  magnitudes.  Hence,  the  latter 
Ms: mb  screening  criterion,  1.25mb-Ms  <  2.20,  is  being  used  as  the  experimental  standard 
criterion  at  the  PIDC,  as  of  1  January  1998. 

2.4  Screening  Probabilities  and  Scores. 

Preliminary  testing  of  the  screening  probabilities  and  scores,  defined  in  Section  1.2.3,  was  also 
performed.  Figure  7  shows  scatterplots  of  depth  versus  (1.25mb-Ms)  estimates  for  all  1996  (upper 
plot)  and  1997  (lower  plot)  events  above  mb  3.5  in  the  REB.  The  circles  depict  events  that  are 
screened  out  at  the  95%  confidence  level,  while  the  triangles  depict  events  that  are  not  screened 
out.  The  95%  confidence  intervals  for  two  events,  including  the  96/06/08  Lop  Nor  explosion,  are 
illustrated  in  the  upper  plot.  Note  that  although  some  of  the  events  have  depth  and/or  Ms:mb 
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estimates  that  are  in  the  “screened  out”  region  of  the  plot,  their  uncertainties  are  large  enough  so 
that  their  95%  confidence  intervals  are  not  entirely  in  the  “screened  out”  region. 

Figure  8  shows  histograms  of  the  screening  probabilities  for  REB  events  above  mb  3.5,  and  with 
unconstrained  depth  and/or  Ms  estimates,  for  1996  (upper)  and  1997  (lower).  The  96/06/08  Lop 
Nor  explosion  has  the  lowest  probability,  by  far,  of  all  the  events  during  these  two  years.  Note  that 
although  a  significant  fraction  of  events  (about  16%)  are  not  screened  out  by  depth  and/or  Ms:mb, 
the  probabilities  indicate  which  events  are  the  most  inconsistent  with  natural  seismic  phenomena. 

As  an  alternative  numerical  representation  of  the  screening  results,  Figure  9  shows  histograms  of 
the  screening  scores  for  the  same  events.  Included  in  the  upper  plot  of  Figure  9  is  the  score  for  the 
96/06/08  Lop  Nor  explosion,  which  is  the  lowest  score  of  all  the  events  in  the  REB  during  1996 
and  1997.  Included  in  the  lower  plot  of  Figure  9  are  the  scores  for  the  93  historical  underground 
nuclear  explosions  at  NTS  and  Semipalatinsk.  Note  that  all  of  these  explosions  at  Lop  Nor,  NTS 
and  Semipalatinsk  have  negative  screening  scores  and,  hence,  none  were  screened  out. 

2.5  Conclusions  and  Recommendations. 

The  experimental  screening  results  presented  here  indicate  that  a  significant  percentage  (about 
77%)  of  the  events  above  mb  3.5  can  be  screened  out  by  the  rather  conservative  criteria,  based  on 
depth,  Ms:mb  and  offshore  location,  and  none  of  the  95  explosions  at  Lop  Nor,  NTS  and 
Semipalatinsk  were  screened  out.  About  20%  of  the  events  above  mb  3.5  had  insufficient  data  with 
which  to  apply  these  screening  criteria.  Three  events  in  the  “insufficient  data”  category  during 
1996  and  1997  were  a  96/07/29  Lop  Nor  explosion,  an  event  in  Southern  India  on  96/1 1/10,  and 
an  event  in  the  Kara  Sea  on  97/08/16.  None  of  these  events  were  screened  out;  however,  these  cases 
illustrate  the  need  for  screening  analyses  based  on  regional  data.  Significant  improvements  in  the 
overall  screening  performance  are  expected  as  the  IMS  network  is  completed  and  once  the 
screening  analysis  based  on  regional  data  is  more  fully  implemented. 

Screening  analysis  based  on  ratios  of  regional  phase  amplitudes  (Pn/Lg  and  Pn/Sn)  have  also  been 
implemented  and  tested.  Applications  to  two  Lop  Nor  explosions  and  other  events  of  interest  are 
described  below  in  Sections  4  and  5.  Fisk  et  al.  (1996b)  describe  other  applications  to  regional 
earthquakes  and  explosions.  Similar  regional  analyses  have  been  tested  for  the  other  IMS  seismic 
stations;  however,  additional  explosion  data  are  needed  before  the  performance  can  be  more  fully 
assessed.  Also,  further  work  is  needed  to  establish  regional  distance  corrections  and  grids  of 
regional  phase  blockage  before  this  analysis  is  ready  to  be  integrated  into  operations  at  the  PIDC. 
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Further  evaluation  of  these  screening  procedures  and  criteria  are  needed,  utilizing  additional 
explosion  data,  in  particular.  To  facilitate  this  evaluation,  work  is  currently  being  conducted,  in 
coordination  with  staff  at  the  Australian  NDC,  to  reprocess  events  in  the  Nuclear  Explosion 
Database  at  the  PIDC.  Details  of  this  database  are  provided  on  the  PIDC  Web  site  (http:// 
www.pidc.org)  under  the  Data  Products  section.  We  plan  to  report  on  this  work  in  the  near  future. 
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Numbers  (upper  plot)  and  percentages  (lower  plot)  of  seismic  events  in  the  REB  per  month  by 
screening  category. 


Percentage  of  Events  Number  of  Events 


Screening  Results  (Numbers)  vs.  mb:  Jul-Dec  97 


Events  Above  mb  Threshold 
Insufficient  Data 
Not  Screened 
Total  Screened  Out 
Screened  by  Depth 
Screened  by  mb-Ms 
Screened  by  Depth  or  mb-Ms 


3.5  4.0  4.5 

Magnitude  Threshold  (mb) 


Screening  Results  (Percentages)  vs.  mb:  Jul-Dec  97 


-  Insufficient  Data 

-  Not  Screened 

-  Total  Screened  Out 

■  -  -  -  Screened  by  Depth 

.  Screened  by  mb-Ms 

-  Screened  by  Depth  or  mb-Ms 


Magnitude  Threshold  (mb) 


Cumulative  numbers  (upper)  and  percentages  (lower)  of  REB  seismic  events  per  screening 
category  as  functions  of  mb  for  July  through  December  1997. 
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Figure  6.  Comparison  of  two  Ms:mb  screening  criteria  for  1996  (upper  plots)  and  1997  (lower  plots).  The 
plots  on  the  left  correspond  to  a  screening  criterion  requiring  that  the  95  %  confidence  interval  for 
mb-Ms  is  less  than  1.2,  while  the  plots  on  the  right  correspond  to  a  screening  criterion  requiring 
that  the  95%  confidence  interval  for  1.25mb-Ms  is  less  than  2.2.  Events  that  are  screened  out  are 
depicted  by  green  circles,  while  those  that  are  not  screened  out  are  depicted  by  red  triangles. 
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Figure  9.  Histogram  plots  of  the  combined  scores  for  all  REB  events  with  mb  above  3.5  and  Ms  and/or 
unconstrained  depth  estimates  during  1996  (upper  plot)  and  1997  (lower  plot).  The  96/06/08  Lop 
Nor  explosion  has  the  lowest  score  of  all  the  events  during  these  two  years. 
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Section  3 

Distance  Corrections  for  Regional  Amplitude  Ratios 


3.1  Introduction. 

The  importance  and  necessity  of  distance  corrections  for  regional  seismic  amplitude  ratios,  in 
order  to  accurately  characterize  regional  seismic  events,  are  thoroughly  discussed  in  previous 
reports  (e.g.,  Bottone  et  al.,  1997;  Fisk  et  ah,  1994,  1995,  1996b). 

Bottone  et  al.  (1997)  described  initial  efforts  to  compute  distance  corrections  for  regional  seismic 
amplitude  ratios,  Pn/Lg  and  Pn/Sn  in  2-4, 4-6  and  6-8  Hz  bands,  for  regions  surrounding  individual 
IMS  seismic  stations.  The  term  “regional”  denotes  events  within  20  degrees  (2220  km)  from  a 
given  station.  In  that  work,  a  simple  power-law  distance  relation  was  used  to  obtain  least-squares 
fits  to  the  regional  data  for  each  IMS  station.  In  that  study,  the  data  used  included  regional  events 
in  the  Reviewed  Event  Bulletin  (REB)  that  had  occurred  between  10  September  1995  and  15 
January  1997.  The  distance  corrections  for  the  individual  stations  were  also  categorized  in  terms 
of  the  number  of  available  regional  events  with  which  to  obtain  the  fits,  as  well  as  whether  the 
distance  corrections  exhibited  anomalous  (or  unexpected)  behavior.  Distance  dependence  with 
either  negative  slopes  or  with  dramatically  steep  slopes  were  considered  anomalous. 

Here  we  present  results  of  continuing  efforts  to  quantify  and  improve  the  distance  corrections  for 
regional  seismic  amplitude  ratios  for  the  existing  IMS  seismic  stations.  In  this  work  we  utilize  an 
improved  parametrization  of  the  distance  dependence,  including  attenuation  and  geometrical 
spreading  terms.  It  will  be  seen  that  the  assumption  of  a  simple  power-law  distance  relation,  used 
in  the  past,  is  too  simplified  to  adequately  fit  the  distance  dependence  over  the  entire  range  of 
regional  distances.  We  also  make  use  of  additional  regional  seismic  data,  collected  at  the  PIDC 
between  10  September  1995  and  1  October  1997.  In  this  work,  we  further  utilize  more  stringent 
restrictions  on  the  quality  of  the  data,  in  terms  of  the  signal-to-noise  ratio  (SNR). 

In  the  remainder  of  this  section,  we  first  provide  descriptions  of  the  regional  seismic  data  that  are 
used.  We  then  describe  the  parameterization  of  the  distance  dependence  and  provide  details  of  the 
worldwide-average  fits  for  all  of  the  regional  data,  collectively.  In  the  following  subsection,  we  use 
the  same  parameterization  to  obtain  region-specific  fits  for  each  individual  IMS  station,  which  are 
compared  to  the  worldwide-average  fits.  We  then  summarize  the  status  of  the  regional  distance 
corrections  on  a  station-by-station  basis.  We  then  describe  and  illustrate  two  approaches  of 
subregionalization,  that  have  been  considered  by  Rodgers  et  al.  (1997)  and  Jenkins  et  al.  (1996),  to 
account  for  variations  of  path  and  subregional  tectonic  effects.  Last,  we  provide  conclusions  and 
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recommendations  based  on  this  work.  An  appendix  to  this  report  provides  further  details  regarding 
the  distance  corrections  for  the  IMS  seismic  stations  with  adequate  regional  data. 

3.2  IMS  Regional  Seismic  Data. 

For  this  analysis,  we  use  all  regional  events  above  mb  3.5  in  the  REB  between  10  September  1995 
and  1  October  1997.  The  restriction  to  events  above  mb  3.5  is  intended  to  avoid  potential 
contamination  from  mining  blasts.  During  that  time  period,  there  were  7,927  regional  events 
recorded  by  the  existing  87  IMS  seismic  stations  (34  primary  and  53  auxiliary)  which  have  at  least 
one  of  the  six  event  characterization  parameters,  Pn/Lg  and  Pn/Sn  in  the  2-4,  4-6,  and  6-8  Hz 
bands,  measured  with  SNR  for  Pn  greater  than  or  equal  to  2.0.  A  more  stringent  SNR  cutoff  of  2.0 
for  both  P  and  S  can  also  be  used;  here  we  follow  Lay  et  al.  (1997). 

It  should  be  noted  that  for  distances  less  than  about  100-300  km  the  regional  Pn  phase  is  not  well 
defined  (Richards,  1997).  Currently  the  PIDC  computes  a  quantity  labeled  as  Pn  by  measuring  the 
maximum  amplitude  in  a  group  velocity  window  which  corresponds  to  Pn  for  distances  greater 
than  about  300  km.  (Section  1.1.5  of  this  report  describes  the  regional  phase  amplitude 
measurements  at  the  PIDC.)  Pmax  and  Smax  amplitude  measurements,  for  distances  up  to  300  km, 
are  planned  to  be  implemented  at  the  PIDC  by  March  1998  to  alleviate  the  problem  (Jepsen,  1998). 
To  avoid  potential  contamination  of  Pn  and  Pg  signals,  events  with  distances  less  than  three 
degrees  (a  conservative  distance  cutoff  recommended  by  Kim,  Lay  and  Walter,  1997)  are  excluded 
from  the  analysis.  Note  that  there  are  only  324  of  the  7,927  events  with  distances  less  than  three 
degrees.  A  sample  of  7,603  regional  events  remains. 

3.3  Worldwide-Averaged  Distance  Corrections. 

Figure  10  contains  semi-log  scatterplots  of  Pn/Lg  and  Pn/Sn  in  the  three  frequency  bands  versus 
distance  from  the  recording  station  for  all  7,603  events.  The  plots  indicate  that  each  P/S  amplitude 
ratio  increases  more  rapidly  as  a  function  of  distance  up  to  about  600  km  and  then  flattens  out  from 
about  600  km  to  2220  km  (20  degrees).  (P/S  is  used  here  to  generically  refer  to  any  of  the  ratios 
Pn/Lg  or  Pn/Sn  in  any  frequency  band.)  To  better  see  the  shape  of  these  curves,  the  data  have  been 
separated  into  34  bins,  each  of  length  0.5  degrees  (55.5  km).  The  mean  of  log(P/S)  in  each  bin  are 
depicted  by  the  square  markers  in  Figure  11,  with  errors  bars  that  are  equal  to  plus  or  minus  the 
estimated  standard  deviation  of  log(P/S)  for  the  data  in  the  bin. 

In  light  of  these  figures,  the  distance  dependence  of  each  P/S  ratio  is  parametrized  as 

log(P/S)  =  a  +  piog(A/A0)  +  yA  ,  (16) 
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where  a,  (3,  and  y  are  parameters  depending  on  the  particular  amplitude  ratio  (Pn/Lg  or  Pn/Sn)  and 
the  frequency  band,/,  A  is  the  epicentral  distance;  and  the  constant  A0  is  a  reference  distance  that 
is  set  at  1500  km,  roughly  the  mean  epicentral  distance  of  the  7,603  regional  events  considered  in 
this  study.  Equation  (16)  is  the  product  of  three  terms:  a  constant,  a  geometrical  spreading  term, 
and  an  attenuation  term.  This  is  a  linear  equation  in  the  unknown  parameters,  a,  (3,  and  y,  which 
can  be  estimated  using  standard  multiple  linear  regression  analysis.  The  regional  amplitude  ratios 
are  distance-corrected  by  applying  the  transformation. 


log(P/S)corrected  =  log  (P/S  )uncorrected  -  a  -  (3log(A/A0)-  yA,  (17) 


where  a ,  (3 ,  and  y  are  the  best  least-squares  estimates  of  a,  (3,  and  y.  Figure  12  contains  plots  of 
Equation  (16)  for  each  amplitude  ratio  using  the  best  least-squares  values  for  a,  [3,  and  y,  based  on 
the  7,603  events.  Superimposed  on  each  plot  are  all  7,603  events  from  Figure  10,  in  gray,  and  the 
binned  means  from  Figure  11.  Also  included  in  Figure  12  is  a  table  giving  the  number  of  regional 
amplitude  measurements,  n,  the  estimated  standard  deviations  of  the  data  both  before  and  after 
distance  corrections,  and  the  best  fit  coefficients,  a,  (3 ,  and  y.  If  we  let  yi  =  log  (P/S  )f, 
i  =  1, ...,  n ,  then  the  estimated  standard  deviation  before  distance  corrections,  s ,  is  defined  by 

52  =  ~i  i>,-502.  (18) 

i  =  1 


where  the  mean,  y ,  is  given  by 


i  =  1 


(19) 


The  estimated  standard  deviation  after  distance  corrections,  sc ,  is  given  by 

n 

sc  =  ~3  X  0>i  -  &  -  p  !°g(  VAo)  -  f  A;)2  •  (20) 

i=  1 


As  can  be  seen  in  the  table,  the  corrected  standard  deviations  are  reduced  by  about  6%  compared 
to  the  uncorrected  values.  If  a  straight-line  fit  is  used,  the  improvement  is  only  about  half  as  much. 
This  estimate  of  the  reduction  in  the  variance  is  not  pronounced  since  there  are  many  more  events 
at  farther  distances  (because  there  is  more  area  at  farther  distances),  where  the  curve  is  relatively 
flat.  Another  measure  of  the  reduction  in  variance  after  distance  corrections,  which  weights  each 
distance  equally,  is  to  compute  the  estimated  standard  deviations  before  and  after  distance 
corrections  for  the  34  binned  values.  In  this  case,  the  average  estimated  standard  deviation  before 
distance  corrections  is  about  0. 18  and,  after  distance  corrections,  about  0.08.  This  is  a  reduction  of 
about  58%,  as  opposed  to  33%  for  the  straight-line  correction.  It  is  not  clear  which  is  the  more 
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Figure  10.  Semi-log  scatterplots  of  regional  amplitude  ratios  versus  distance  for  7,603  events  above  mb  3.5 
with  SNR  for  Pn  greater  than  or  equal  to  2.0. 
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Figure  11.  Semi-log  plots  of  mean  (log)  regional  amplitude  ratios  versus  distance  for  7,603  events.  Data  is 
separated  into  bins  of  length  55.5  km  (0.5  degrees).  Error  bars  are  equal  to  plus  or  minus  one 
standard  deviation. 
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appropriate  way  to  estimate  the  standard  deviations,  so  below  the  more  conservative  computation 
which  treats  each  event  with  equal  weight  is  used. 

Figure  12  indicates  that  a  straight  line  will  not  approximate  the  distance  dependence  well  for  all 
distances,  due  to  the  steeper  slopes  at  distances  less  than  about  600  km.  Also,  the  corrections  will 
be  most  important  for  distances  less  than  about  600  km.  The  flatness  of  the  curves  beyond  about 
600  km  indicates  that  distance  corrections  are  not  as  necessary  for  events  at  these  distances. 

3.4  Region-Specific  Distance  Corrections. 

Equation  (16)  provides  good  fits  to  the  combined  regional  data  from  all  IMS  stations,  although 
there  is  considerable  scatter  in  the  data  due,  in  large  part,  to  regional  variations.  Equation  (16)  is 
now  used  to  fit  the  regional  data  from  individual  IMS  stations.  There  are  currently  50  IMS  stations 
that  have  30  or  more  regional  events  (with  SNR  >  2.0 )  with  which  to  estimate  the  fits. 

As  an  example,  plots  of  the  best  least-squares  fits  (solid  curves)  for  1 12  regional  events  recorded 
by  station  NIL  are  shown  in  Figure  13.  Also  shown  in  the  plots,  for  comparison,  are  the  worldwide 
average  fits  (dashed)  from  Figure  12.  For  reference,  the  plots  also  include  the  P/S  values  for  two 
underground  explosions  at  the  Lop  Nor  test  site  during  1996.  These  two  events  are  enclosed  in  a 
circle  and  are  outliers  from  the  rest  of  the  sample  population.  For  a  full  account  of  the 
characterization  of  the  Lop  Nor  explosions  see  Fisk  et  al.  (1996a).  The  table  included  in  Figure  13 
gives  the  estimated  standard  deviation  of  the  training  set  both  before,  s,  and  after,  sc,  distance 
corrections,  in  addition  to  the  best-fit  parameters  for  the  distance-correction  curves.  The  average 
reduction  in  standard  deviation  is  about  10%.  It  is  clear  from  the  plots  that  the  distance-correction 
curves  for  NIL  are  nearly  identical  to  the  worldwide  average. 

Two  examples  of  stations  whose  distance  corrections  significantly  reduce  the  variance  are  PDAR 
(Figure  14),  with  93  regional  events,  and  PDY  (Figure  15),  with  59  regional  events.  The  average 
reduction  in  estimated  standard  deviation  is  24%  for  PDAR  and  53%  for  PDY.  The  curves  for  PDY 
are  quite  different  than  the  worldwide-average  curves,  especially  at  short  distances.  This  example 
shows  the  importance  of  determining  appropriate  corrections  for  each  region  individually.  Before 
distance  corrections,  the  average  standard  deviations  are  greater  than  0.5,  while  after  distance 
corrections,  the  standard  deviations  are  mostly  less  than  0.3,  a  relatively  small  value,  thus 
improving  the  event  characterization  capability  in  this  region. 
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Figure  12.  Best  least-square  fit  curves  for  Pn/Lg  and  Pn/Sn  in  the  2-4,  4-6  and  6-8  Hz  bands,  based  on  7,603 
regional  events.  Superimposed  are  the  scatterplots  as  in  Figure  10  and  the  binned  means  as  in 
Figure  11.  The  table  gives  the  number  of  measurements  for  each  amplitude  ratio,  the  estimated 
standard  deviation  of  the  data  before  and  after  distance  corrections,  and  the  best-fit  coefficients. 
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Figure  13.  Best  least-square  fits  (solid  curves)  for  112  regie 
are  the  worldwide  averages  as  in  Figure  12.  The 
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Figure  15.  Best  least-square  fits  (solid  curves)  for  59  regional  events  recorded  by  PDY.  The  dashed  curves  are 
the  worldwide  averages  as  in  Figure  12. 
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Figure  16.  Best  least-square  fits  (solid  curves)  for  211  regional  events  recorded  by  TXAR.  The  dashed  curves 
are  the  worldwide  averages  as  in  Figure  12. 
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As  another  example,  Figure  16  shows  the  data  and  distance-correction  curves  for  TXAR  (with  211 
regional  events).  The  curves  are  quite  similar  to  the  worldwide-average  curves  for  distances 
beyond  about  800  km,  but  differ  for  closer  distances.  There  are,  however,  no  events  for  TXAR  at 
distances  less  than  about  600  km  which  satisfy  the  event  selection  criteria.  Such  data  are  needed  in 
order  to  improve  the  fits  for  distances  less  than  800  km.  The  estimated  distance  corrections  do, 
however,  provide  good  fits  to  the  available  data  and  the  scatter  in  the  data  is  relatively  small. 

Figure  17  shows  the  amplitude  ratios  and  distance-correction  curves  for  1,266  events  recorded  by 
WRA,  all  of  which  have  distances  greater  than  1000  km.  Thus,  there  are  no  data,  satisfying  the  data 
quality  restrictions,  that  can  be  used  to  quantify  the  distance  dependence  at  close  distances.  Note 
that  for  any  event  at  the  far  regional  distances,  where  there  is  much  training  data,  the  distance 
dependence  is  relatively  insignificant  compared  to  the  scatter  in  the  data.  Additional  data  at  closer 
distances  are  needed  to  be  able  to  characterize  events  in  that  distance  range. 

A  station  with  unusual  distance  dependence  is  LPAZ,  with  282  regional  events  (Figure  18).  Each 
of  the  best  fit  distance-correction  curves  differ  significantly  from  the  respective  worldwide 
averages  and  there  is  significant  scatter  in  the  data.  The  complicated  subregional  effects  near  LPAZ 
are  described  below  to  shed  some  light  on  the  underlying  reasons  for  this  behavior. 

3.5  Status. 

The  status  of  the  regional  distance  corrections  is  now  summarized  for  each  existing  IMS  seismic 
station.  For  the  purposes  of  event  screening,  the  important  quantities  are  the  standard  deviations  of 
the  corrected  data,  relative  to  the  separation  of  the  means  for  explosions  and  earthquakes.  It  is  also 
important  to  have  sufficient  high-quality  data  over  a  sufficient  distance  range  with  which  to  obtain 
the  fits.  As  illustrated  above,  these  considerations  are  often  most  important  for  events  at  distances 
less  than  about  600  km  where  the  curves  typically  exhibit  the  strongest  dependence  on  distance. 

Assessing  the  utility  of  the  distance  corrections  would  be  easier  if  regional  data  existed  for  both 
earthquakes  and  explosions  at  the  same  station.  Unfortunately,  explosion  data  exist  for  relatively 
few  stations.  Thus,  our  objective  is  to  reduce  the  variance  as  much  as  possible,  using  distance 
corrections  and  other  procedures,  such  as  dividing  the  data  into  subsets  of  reduced  variance  (see 
Lay  et  al.,  1997,  and  below).  Stations  that  do  have  regional  explosion  data  can  then  be  used  as  a 
guide  in  assessing  how  small  the  variance  should  be  in  order  to  be  useful  in  event  screening. 
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Figure  17.  Best  least-square  fits  (solid  curves)  for  1266  regional  events  recorded  by  WRA.  The  dashed  curves 
are  the  worldwide  averages  as  in  Figure  12. 
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Figure  18.  Best  least-square  fits  (solid  curves)  for  282  regional  events  recorded  by  LPAZ.  The  dashed  curves 
are  the  worldwide  averages  as  in  Figure  12. 


Here,  we  characterize  the  training  sets  in  two  ways.  First,  since  the  distance-correction  curves 
given  in  Equation  (16)  generally  work  quite  well,  given  sufficient  data,  we  categorize  stations 
simply  by  the  number  of  events  with  at  least  one  regional  amplitude  ratio  with  SNR  >  2.0 .  Table 
5  lists  each  of  the  34  IMS  primary  stations  and  their  corresponding  category.  A  similar  list  appears 
in  Table  6  for  the  53  IMS  auxiliary  stations.  Stations  listed  as  Category  1  have  at  least  30  regional 
seismic  events  with  SNR  >  2.0,  while  those  listed  as  Category  2  do  not.  (The  station  names  are 
those  currently  used  at  the  PIDC.)  Also  included  in  each  table  are  the  number  of  regional  events 
recorded  at  each  station  between  10  Sep  1995  and  1  Oct  1997  with  SNR  >2.0.  The  column 
labeled  smin  gives  the  minimum  of  the  six  distance-corrected  standard  deviations.  The  column 
labeled  dmin  gives  the  minimum  distance  from  event  to  station  of  all  events  in  the  training  set. 

Summarizing  Table  5  and  Table  6,  there  are  50  IMS  stations  (27  primary  and  23  auxiliary)  that 
currently  have  adequate  data  with  which  to  estimate  region-specific  distance  corrections  and  37 
stations  (7  primary  and  30  auxiliary)  that  do  not.  This  situation  will  improve  significantly  over 
time,  as  more  regional  data  are  acquired  from  the  IMS  seismic  stations. 

As  a  second  method  of  categorization,  we  determine  the  number  of  stations  for  which  the  standard 
deviation  after  distance  corrections,  sc,  is  below  a  specified  level,  s0,  for  all  six  regional  amplitude 
ratios.  As  discussed  above,  the  appropriate  value  for  s0  can  be  determined  precisely  only  if  the 
variance  and  difference  in  the  means  of  both  earthquakes  and  explosions  are  known,  and  is  most 
likely  highly  station  dependent.  Without  this  information,  we  simply  chart  the  number  of  stations 
for  various  values  of  Sq,  summarized  in  Table  7.  As  can  be  seen  in  the  table,  over  80%  of  the  50 
stations  with  adequate  data  have  standard  deviations  after  distance  corrections  less  than  0.5,  for  all 
six  regional  amplitude  ratios,  and  over  95%  less  than  0.6. 

Figure  19  shows  histograms  of  the  number  of  stations  with  standard  deviations  between  specified 
values,  both  before  and  after  distance  corrections.  Figure  19  shows,  for  example,  that  the  number 
of  stations  which  have  standard  deviations  less  than  0.5  for  all  six  regional  amplitude  ratios  before 
distance  corrections  is  only  27  (54%)  as  opposed  to  41  (82%)  after  distance  corrections.  Only  12 
stations,  as  opposed  to  25,  have  all  six  standard  deviations  less  than  0.4,  and  only  42  stations,  as 
opposed  to  48,  have  all  six  standard  deviations  less  than  0.6  before  distance  corrections. 
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Table  5.  Status  of  regional  distance  corrections  for  IMS  primary  seismic  stations. 


Table  6.  Status  of  regional  distance  corrections  for  IMS  primary  seismic  stations. 


Station 

Category 

#  Events 
(SNR  >  2.0) 

smin 

dmin 

AAb 

‘  "  2 

N/A 

356 

AF1 

i 

71 

0.385 

336 

ALQ 

i 

61 

0.162 

852 

AQU 

2 

3 

N/A 

505 

ARU 

i 

30 

0.224 

1712 

BBB 

2 

27 

0.352 

333 

BORG 

2 

9 

0.213 

383 

CTA 

i 

215 

0.243 

986 

DAV 

2 

20 

0.177 

404 

DAVOS 

i 

179 

0.268 

351 

DLBC 

i 

104 

0.272 

355 

EIL 

2 

27 

0.264 

573 

ERA 

1 

158 

0.252 

645 

ELK 

1 

44 

0.371 

379 

FITZ 

2 

25 

0.196 

1310 

GNI 

2 

1 

N/A 

343 

HFS 

i 

111 

0.242 

569 

HNR 

i 

49 

0.332 

357 

INK 

i 

118 

0.358 

355 

JCJ 

2 

0 

N/A 

N/A 

JHJ 

2 

0 

N/A 

N/A 

JKA 

2 

0 

N/A 

N/A 

JNU 

2 

■■■■ 

N/A 

N/A 

JOW 

2 

1 

N/A 

1418 

JTS 

2 

0 

N/A 

N/A 

KDAK 

2 

10 

0.322 

1310 

KIEV 

i 

46 

0.336 

438 

KVAR 

i 

36 

0.198 

335 

LSZ 

2 

0 

N/A 

N/A 

MLR 

2 

0 

N/A 

N/A 

MSEY 

2 

0 

N/A 

N/A 

NEW 

i 

50 

0.230 

362 

NIL 

i 

112 

0.245 

394 

NNA 

i 

37 

0.443 

341 

NWAO 

2 

18 

0.370 

1316 

OBN 

1 

36 

0.225 

381 

PARD 

2 

9 

0.163 

370 

PFO 

1 

53 

0.280 

444 

PMG 

1 

96 

0.219 

340 

PTGA 

2 

0 

N/A 

N/A 

RAR 

2 

28 

0.207 

1367 

RPN 

2 

3 

N/A 

597 

SADO 

2 

4 

0.089 

929 

SDV 

2 

0 

N/A 

N/A 

SFJ 

2 

11 

0.373 

833 

SNZO 

1 

71 

0.266 

376 

SPITS 

1 

77 

0.346 

384 

SUR 

2 

16 

0.209 

711 

TKL 

2 

7 

0.241 

1889 

TSUM 

2 

7 

0.054 

878 

ULN 

i 

33 

0.150 

465 

VRAC 

i 

58 

0.206 

414 

YAK 

2 

0 

N/A 

N/A 

Table  7.  Number  of  stations  with  distance-corrected  standard  deviations  less  than  s0 


So 

Primary 

r>(sc<s0) 

Auxiliary 

n(sc<s0) 

Total 

n(sc<so) 

0.1 

0 

0 

0 

0.2 

0 

0 

0 

0.3 

3 

3 

0.4 

16 

9 

25 

0.5 

23 

18 

41 

0.6 

26 

22 

48 

0.7 

26 

23 

49 

0.8 

27 

23 

50 

0.9 

27 

23 

50 

1.0 

27 

23 

50 

3.6  Subregional  Improvements. 

The  examples  given  above  for  TXAR,  WRA  and  LPAZ  show  that  distance  corrections  alone  do  not 
always  appreciably  reduce  the  variance  of  the  training  data,  due  to  subregional  or  path-specific 
variations.  In  some  cases,  subsets  of  a  training  set,  i.e.,  data  restricted  to  a  particular  subregion,  can 
provide  better  event  characterization.  Regional  amplitude  ratios  of  events  in  such  a  subregion 
could  have  a  standard  deviation  less  than  that  of  the  entire  regional  set.  This  reduced  variance 
could  depend  on  the  properties  of  the  locations  of  the  events  or  on  the  path  differences  from  the 
events  to  the  station,  or  both.  Such  variations  should  be  treated  in  order  to  provide  more  effective 
regional  event  screening.  Tectonic  subregion  and  azimuthal  dependence  are  examined  here  for 
TXAR,  WRA  and  LPAZ. 

3.6.1  Treatment  of  Subregional  Tectonic  Variations. 

Figures  20-22  show  the  locations  of  the  regional  events  recorded  by  LPAZ,  TXAR  and  WRA, 
respectively,  on  maps  of  each  region.  Each  map  shows  the  location  of  the  station,  as  well  as  a 
tectonic  grid,  with  two-degree  resolution,  that  was  established  by  Oli  Guudmundsson  of  the 
Australian  National  University,  and  is  described  further  by  Jenkins  et  al.  (1996).  Each  grid  is 
divided  into  up  to  eight  distinct  tectonic  subregions,  color-coded  in  the  following  way: 
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ALL:  Pn/Lg(2-4  Hz)  ALL:  Pn/Sn(2-4  Hz) 


Figure  19.  Histograms  showing  the  number  of  stations 
distance  corrections. 
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cyan: 

Young  ocean  (<25  Myrs) 

blue: 

Intermediate  ocean  (25-100  Myrs) 

black: 

Old  ocean  (>100  Myrs) 

green: 

Cenozoic  (alpine) 

magenta: 

Tectonic 

red: 

Paleozoic  (ca  150-800  Myrs) 

yellow: 

Proterozoic  (800-1700  Myrs) 

white: 

Archean  (>1700  Myrs) 

These  figures  depict  the  distribution  of  events  relative  to  tectonic  subregion  and  station  location. 
For  example,  most  regional  events  recorded  by  LPAZ  occur  along  the  Andes,  either  to  the 
northwest  of  the  station  or  to  the  south,  either  in  the  tectonic  (magenta),  Cenozoic  (green),  or 
intermediate  ocean  (blue)  subregions.  The  regional  amplitude  ratios  for  LPAZ  are  plotted  versus 
distance  in  Figure  23  for  each  subregion  on  the  same  plot,  along  with  the  fits  of  the  distance 
dependence,  color-coded  as  above.  Events  from  each  subregion  do  have  different  distance 
dependence,  with  only  the  tectonic  (magenta)  events  consistent  with  the  worldwide  averages.  The 
intermediate  ocean  (blue)  events  show  typical  behavior  in  the  2-4  Hz  frequency  band,  but  a 
decreasing  dependence  with  distance  at  the  higher  frequencies.  The  events  in  the  Cenozoic  (green) 
subregion  have  a  relatively  flat  distance  dependence.  The  standard  deviations  of  the  distance- 
corrected  data  for  Pn/Sn(6-8  Hz),  for  example,  is  0.30  for  the  tectonic,  0.30  for  the  Cenozoic,  and 
0.27  for  the  intermediate  ocean  subregions,  as  compared  to  0.36  for  all  the  data.  Thus,  in  this  case, 
there  is  something  to  be  gained  by  making  this  subregional  division.  However,  additional  data  and 
further  evaluation  are  needed  for  LPAZ  to  validate/improve  these  corrections. 

Similar  plots  are  shown  in  Figure  24  for  TXAR  for  four  tectonic  subregions:  tectonic  (magenta), 
Paleozoic  (red),  young  ocean  (cyan),  and  intermediate  ocean  (blue).  The  individual  distance- 
correction  curves  all  exhibit  similar  behavior  at  far  regional  distances;  the  erratic  behavior  at 
distances  less  than  500-1000  km  is  due  to  a  lack  of  data  in  this  distance  range.  Unlike  LPAZ,  there 
is  hardly  any  reduction  in  standard  deviation  for  data  confined  to  any  of  the  tectonic  subregions, 
except  for  the  intermediate  ocean  events.  However,  as  was  pointed  out  earlier,  the  variance  of  all 
the  TXAR  regional  events  is  already  relatively  small.  Azimuthal  variations  in  the  regional  data 
recorded  by  TXAR  will  be  studied  further  below. 
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Figure  20.  Regional  event  locations  and  tectonic  grid  for  LPAZ. 


Longitude 


Figure  21.  Regional  event  locations  and  tectonic  grid  for  TXAR. 
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Figure  22.  Regional  event  locations  and  tectonic  grid  for  WRA. 

3.6.2  Treatment  of  Azimuthal  Variations. 

Another  approach  involves  dividing  regions  into  azimuthal  sectors.  Successful  application  of  this 
approach  for  ABKT  is  reported  by  Rodgers  et  al.  (1997).  Events  measured  at  similar  distances  but 
different  azimuths  provide  information  about  path  differences.  In  Figures  25-27,  the  azimuthal 
variations  for  the  six  P/S  ratios  are  plotted  for  LPAZ,  TXAR  and  WRA,  respectively.  Azimuths  are 
measured  clockwise  from  north.  The  data  for  LPAZ  and  TXAR  have  been  divided  into  azimuthal 
bins  of  size  18  degrees  and,  for  WRA,  into  bins  of  5  degrees.  The  solid  line  in  each  plot  represents 
the  estimated  mean  in  each  bin,  and  the  errors  bars  are  plus  and  minus  one  estimated  standard 
deviation  of  the  data  in  the  bin.  The  two  horizontal  dashed  lines  in  each  plot  correspond  to  plus  and 
minus  one  standard  deviation  relative  to  the  mean  of  the  entire  data  set  for  that  station. 

For  LPAZ,  the  large  error  bars  for  azimuths  between  30  and  120  degrees  are  most  likely  due  to  the 
lack  of  data  at  these  angles  (see  Figure  20).  Figure  20  also  shows  that  most  events  occur  to  the 
south  (150-210  degrees)  and  northwest  (270-330  degrees)  relative  to  LPAZ.  For  most  of  the  P/S 
ratios,  the  means  for  northwestern  azimuths  are  somewhat  smaller  than  for  southern  azimuths,  as 
are  the  standard  deviations.  There  are  two  azimuthal  bins,  at  about  135  degrees  and  260  degrees, 
which  have  significantly  smaller  standard  deviations,  as  compared  to  the  dashed  lines.  However, 
not  every  azimuthal  range  has  standard  deviations  less  than  the  overall  standard  deviation. 
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Figure  23. 


Figure  24. 


LPAZ:  Pn/Lg(2-4  Hz) 


LP.4Z:  Pn/Lg(4-6  Hz) 


LPAZ:  Pn/Lg(6-8  Hz) 


LPAZ:  Pn/3n(2-  4  Ilz) 


LPAZ:  Pn/Sn(4-  6  Hz) 


LPAZ:  Pn/3n(6-8  Hz) 
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Best  least-square  fit  curves  for  LPAZ;  magenta  for  events  in  tectonic  subregions,  green  for 
Cenozoic,  and  blue  for  intermediate  ocean. 


TXAR:  Pn/Lg(2  — 4  Hz)  TXAR:  Pn/Sn(2-4  Hz) 


TXAR:  Pn/Lg(6 -8  Hz) 


TXAR:  Pn/Sn(6-8  Hz) 


Best  least-square  fit  curves  for  TXAR;  magenta  for  events  in  tectonic  subregions,  red  for 
Paleozoic,  cyan  for  young  ocean,  and  blue  for  intermediate  ocean. 
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Figure  25.  Azimuthal  variations  at  LPAZ.  Error  bars  are  plus  or  minus  one  standard  deviation. 


Figure  26  shows  that  the  azimuthal  variations  of  the  standard  deviations  for  TXAR  do  not  vary 
appreciably  compared  to  all  of  the  data.  The  means  of  the  azimuthal  bins,  however,  do  exhibit  a 
significant  variation  near  240  degrees,  for  events  that  are  mostly  in  the  Gulf  of  California, 
especially  at  lower  frequencies.  Treatment  of  the  azimuthal  variation  can  be  used  to  further  reduce 
the  variance  of  the  regional  amplitude  ratios  for  TXAR. 

WRA  appears  to  be  a  good  station  to  study  azimuthal  variations  because  there  are  numerous  events 
with  similar  distances  from  the  station,  with  azimuths  ranging  from  approximately  -60  to  about  60 
degrees,  as  shown  in  Figure  22.  Figure  27  shows  azimuthal  variations  in  5  degree  bins.  The  means 
are  quite  uniform  over  the  entire  120  degree  range.  There  is  some  variation  in  standard  deviation, 
the  smallest  standard  deviation  having  about  half  the  value  of  the  largest  standard  deviation.  It  is 
not  the  case,  as  would  be  most  useful,  that  all  standard  deviations  in  the  azimuthal  bins  are  less 
than  the  values  for  all  the  WRA  data. 
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Figure  26. 


Figure  27. 


TXAR:  Pn/Lg(2-4  Hz) 


TXAR:  Pn/Sn(2-4  Hz) 


TXAR:  Pn/Lg(4-6  Hz) 


TXAR:  Pn/Lg(6-8  Hz) 


TXAR:  Pn/Sn(4-6  Hz) 


TXAR:  Pn/Sn(6-8  Hz) 
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Azimuthal  variations  at  TXAR.  Error  bars  are  plus  or  minus  one  standard  deviation. 


WRA:  Pn/Lg(2-4  Hz) 


WRA:  Pn/Sn(2-4  Hz) 


WRA:  Pn/Lg(4-6  Hz) 


WRA:  Pn/Lg(6-8  Hz) 


WRA:  Pn/Sn(4-6  Hz) 


WRA:  Pn/Sn(6-8  Hz) 


Azimuthal  variations  at  WRA.  Error  bars  are  plus  or  minus  one  standard  deviation. 
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3.7  Conclusions  and  Recommendations. 

In  this  section,  we  described  ongoing  work  to  improve  distance  corrections  for  regional  phase 
amplitude  ratios,  utilizing  7,603  regional  events  recorded  by  the  IMS  seismic  network.  The 
distance  corrections  include  both  geometrical  spreading  and  attenuation  terms,  and  provide  a  much 
better  approximation  of  the  distance  dependence  than  does  a  simple  power-law  parametrization. 
Roughly  50  of  the  existing  87  IMS  seismic  stations  have  adequate  data  (i.e.,  at  least  30  events  with 
SNR  >  2.0 )  with  which  to  compute  the  distance  corrections.  Explicit  details  of  the  distance 
corrections  for  these  stations  are  provided  in  the  appendix  to  this  report. 

Anomalous  distance  dependence  at  some  stations  can  be  attributed,  in  part,  to  variations  in 
propagation  paths  and  properties  of  the  near-source  media.  Treatments  of  azimuthal  variations  and 
tectonic  subregions  were  illustrated  for  LPAZ,  TXAR  and  WRA,  showing  possible  improvements 
that  can  be  made  by  dividing  training  sets  into  subregions,  each  with  a  smaller  variance  than  the 
entire  set.  Another  promising  approach  by  Lay  et  al.  (1997),  treating  other  geophysical  properties 
along  the  path,  such  as  mean  topography,  crustal  thickness,  sediment  thickness,  and  roughness,  has 
been  applied  to  station  ABKT.  These  and  other  approaches  need  to  be  examined  further. 

Another  approach  that  could  alleviate  the  need  for  distance  corrections  as  IMS  regional  seismic 
data  accumulate  is  to  compare  an  event  in  question  to  other  events  that  are  close  in  distance.  The 
definition  of  “close”  could  mean  having  a  similar  distance  from  the  station,  similar  azimuth  and/or 
tectonic  subregion.  A  better  definition  of  “close”  would  require  all  training  events  to  be  within  a 
certain  distance  of  an  event  being  tested.  This  would  require  much  more  data  so  as  to  have  enough 
training  events  at  roughly  the  same  location.  In  either  case,  it  is  a  non-trivial  problem  to  quantify 
the  notion  of  close.  Basically,  it  must  be  determined  if  it  is  better  to  include  a  potential  training 
event  at  a  given  distance  away  from  the  test  event  than  to  not  include  it  in  the  set.  This  depends  on 
how  different  the  distribution  of  events  are  at  further  locations.  If  this  were  known,  then  a 
correction  could  be  made;  however,  in  these  cases,  the  dependence  on  location  is  not  known  and  it 
is  often  difficult  to  approximate.  Thus,  the  problems  of  determining  the  dependence  of  regional 
amplitude  ratio  distributions  on  location  and  what  is  meant  by  close  are  intimately  related. 

As  more  IMS  data  become  available  at  the  prototype  IDC,  we  plan  to  progressively  improve  and 
validate  the  regional  distance  corrections.  We  are  also  developing  grids  of  blockage  and  strong 
attenuation  for  regional  seismic  phases.  We  plan  to  further  develop,  test  and  evaluate  these  grids. 
We  recommend  that  each  IMS  station  be  studied  individually  to  determine  the  best  approach  of 
dividing  a  particular  regional  training  set  into  subsets  of  minimized  variance.  We  also  recommend 
that  ground-truth  data  for  historical  explosions  be  processed  according  to  standard  PIDC 
techniques  and  used  in  the  evaluation  of  the  regional  distance  corrections  and  screening  analysis. 
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Section  4 

Analysis  of  the  10  November  1996  Southern  India  Event 

4.1  Introduction. 

Based  on  automated  and  analyst-reviewed  processing  at  the  PIDC,  an  mb  4.17  seismic  event 
occurred  on  10  November  1996  (9611 10)  at  approximately  09:00:00.4  UTC,  18.087°  latitude  and 
76.506°  longitude  in  Southern  India.  The  event  was  approximately  35.9°  from  ZAL,  21.2°  from 
CMAR,  and  15.8°  from  NIL.  Signals  were  also  detected  by  19  other  stations.  The  focal  depth  of 
the  event  was  constrained  to  the  surface  in  the  reviewed  location  analysis  and,  hence,  it  cannot  be 
assumed  that  it  was  a  deep  natural  event.  Furthermore,  there  were  no  long-period  surface  waves 
detected.  Thus,  screening  criteria  based  on  depth  and  mb-Ms  could  not  be  applied. 

Using  seismic  waveforms  from  station  NIL  in  Pakistan,  we  compare  this  event  to  previous  events 
in  the  region,  including  presumed  earthquakes  and  two  underground  nuclear  explosions  at  the  Lop 
Nor  test  site  in  China.  (Fisk  et  al.,  1996a,  provide  a  more  detailed  discussion  of  the  Lop  Nor 
events.)  Since  the  961 1 10  event  was  not  necessarily  deep  and  did  not  have  an  estimate  of  Ms,  we 
use  high-frequency  regional  measurements  of  Pn/Lg  and  Pn/Sn  to  make  the  comparisons.  We 
apply  a  statistical  outlier  test  to  make  a  multivariate  comparison  to  previous  events,  as  described  in 
Section  4.3.  Last,  we  provide  some  conclusions  regarding  this  event. 

4.2  Data. 

Figure  28  shows  the  locations  of  station  NIL,  the  961 1 10  event,  two  known  underground  nuclear 
explosions  at  the  Lop  Nor  test  site  in  China,  and  other  regional  seismic  events  (within  20  degrees 
from  NIL)  that  were  considered  in  the  analysis.  These  other  regional  events  are  all  mb  3.5  or  above, 
and  are  presumed  to  be  earthquakes.  Their  paths  to  NIL  represent  a  relatively  broad  range  of 
epicentral  distances  and  azimuths,  including  another  event  in  India  with  a  similar  path  to  NIL. 

Figure  29  shows  raw  and  bandpass-filtered  seismograms  on  the  sz  channel  for  the  96/1 1/10  event 
at  NIL.  Preliminary  inspection  suggests  that,  although  the  signal-to-noise  ratio  was  relatively  low, 
the  961110  event  generated  comparable  high-frequency  Pn,  Lg  and  Sn  phase  amplitudes.  Pn,  Lg 
and  Sn  maximum  amplitude  values  were  computed  at  the  PIDC.  Table  8  summarizes  the  Pn/Lg 
and  Pn/Sn  values  in  the  2-4, 4-6  and  6-8  Hz  bands,  before  distance  corrections  were  applied. 

For  comparison,  Figure  30  shows  the  bandpass-filtered  (4-8  Hz)  seismograms  at  station  NIL  for 
the  961110  event  and  the  960608  (mb  5.69)  and  960729  (mb  4.71)  Lop  Nor  nuclear  explosions. 
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Note  that  the  Lop  Nor  explosions  both  occurred  approximately  14.4°  from  station  NIL,  roughly  the 
same  distance  (but  different  azimuth)  as  the  961110  event  in  India. 


Figure  28.  Locations  of  Auxiliary  station  NIL  and  seismic  events  used  in  the  event  characterization  analysis 
of  the  961110  Southern  India  event. 


Table  8.  Pn/Lg  and  Pn/Sn  values  for  the  961110  event  based  on  waveforms  recorded  by  NIL. 


2-4  Hz 

4-6  Hz 

6-8  Hz 

Pn/Lg 

1.20 

0.87 

0.65 

Pn/Sn 

1.30 

1.66 

1.09 
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Figure  30.  Bandpass  filtered  (4-8  Hz)  seismograms  recorded  by  station  NIL  for  the  960608  and  960729  Lop 
Nor  explosions  and  the  961110  India  event  (bottom  trace).  By  comparison,  the  961110  India  event 
has  considerably  lower  high-frequency  P-to-S  energy  than  the  Lop  Nor  explosions.  Its  Pn  phase  is 
also  far  less  impulsive  as  compared  to  the  Lop  Nor  explosions. 


Note  that  both  Lop  Nor  explosions  had  depth  estimates  constrained  to  the  surface.  The  960608  Lop 
Nor  event  has  a  99%  confidence  interval  for  mb-Ms  entirely  above  a  threshold  of  1 .2,  while  no  Ms 
measurements  were  available  for  the  960729  Lop  Nor  explosion.  Thus,  as  for  the  961110  India 
event,  regional  analysis  was  required  in  order  to  characterize  the  960729  Lop  Nor  event. 

Figure  31  shows  the  Pn/Lg  (left)  and  Pn/Sn  (right)  values  versus  distance  in  the  various  frequency 
bands  for  the  961 1 10  Southern  India  event,  the  two  Lop  Nor  explosions  and  the  other  events  in  the 
same  region,  which  are  presumed  to  be  earthquakes.  Figure  32  shows  plots  of  the  distance- 
corrected  Pn/Lg  and  Pn/Sn  values,  computed  for  the  NIL  recordings  in  the  2-4,  4-6  and  6-8  Hz 
bands.  The  values  for  the  961110  Southern  India  event  appear  to  be  consistent  with  other  regional 
seismic  activity,  and  considerably  different  that  the  Lop  Nor  nuclear  explosions. 


NIL:  Pn/L<g(2-4  Hz) 


Figure  31.  Regional  Pn/Lg  and  Pn/Sn  values  versus  distance,  computed  for  station  NIL  in  the  2-4,  4-6  and  6- 
8  Hz  bands,  for  presumed  earthquakes,  two  Lop  Nor  underground  nuclear  explosions  and  the 
961110  Southern  India  event.  The  lines  depict  the  empirical  distance  dependence  that  is  used  to 
correct  for  regional  attenuation. 
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4.3  Regional  Population  Analysis. 

To  compare  the  961110  event  to  previous  seismic  activity  observed  by  NIL,  we  performed  a 
multivariate  outlier  test,  based  on  the  likelihood  ratio,  in  which  the  961110  event  was  tested  as  an 
outlier  of  the  regional  earthquake  group.  Pn/Lg  and  Pn/Sn  in  the  2-4,  4-6  and  6-8  Hz  bands  were 
used  as  event  characterization  parameters.  For  comparison,  we  also  tested  each  of  the  two  Lop  Nor 
underground  nuclear  explosions.  Using  the  leave-one-out  procedure,  we  further  tested  each  of  the 
earthquakes  as  outliers  of  the  remaining  earthquake  group.  (Fisk  et  al.,  1993,  1994,  1995,  1996a, 
1996b,  describe  the  outlier  procedure  and  previous  applications  in  greater  detail.) 

Figure  33  shows  the  results  of  this  analysis.  The  distribution  shown  corresponds  to  that  of  the 
likelihood  ratio  obtained  from  the  multivariate  discriminant  values  for  the  training  set.  Each 
marker  corresponds  to  the  P-values  for  the  individual  events  being  tested.  The  upper  legend 
associates  the  various  markers  with  the  different  types  of  events.  The  P-value  is  the  probability  that 
a  particular  event  is  called  an  outlier,  given  its  discriminant  values,  if  it  in  fact  belongs  to  the 
training  group.  That  is,  it  is  the  area  under  the  distribution  of  the  likelihood  ratio  to  the  left  of  the 
likelihood  ratio  value  for  the  particular  event.  A  small  P-value  indicates  that  the  event  is 
inconsistent  with  the  training  set.  The  vertical  line  shown  corresponds  to  the  threshold  of  the 
outlier  test  for  a  0.01  significance  level.  Events  whose  P-values  fall  to  the  left  of  this  line  are 
considered  outliers  at  the  0.01  significance  level. 

Note  that  the  P-values  for  the  Lop  Nor  nuclear  explosions  are  both  less  than  10'10.  This  indicates 
that  the  regional  Pn/Lg  and  Pn/Sn  values  for  the  Lop  Nor  explosions  are  extremely  inconsistent 
with  those  for  the  regional  earthquake  group.  The  P-value  for  the  961110  event  is  0.60.  Thus,  it  is 
not  considered  to  be  an  outlier  of  the  regional  population. 

4.4  Conclusions. 

This  investigation  of  the  961110  Southern  India  event  indicates  that  its  regional  high-frequency 
Pn/Lg  and  Pn/Sn  characteristics  are  consistent  with  regional  earthquakes  seen  by  station  NIL. 
Furthermore,  its  regional  Pn/Lg  and  Pn/Sn  characteristics  are  considerably  different  than  those  of 
two  Lop  Nor  underground  nuclear  explosions.  The  regional  population  analysis  indicates  that  there 
is  insufficient  seismic  evidence  to  reject  this  event  as  belonging  to  the  earthquake  group  in  this 
region,  based  on  the  regional  amplitude  ratios  and  training  events  used.  By  comparison,  the 
regional  population  analysis  clearly  indicates  that  both  of  the  Lop  Nor  underground  nuclear 
explosions  are  outliers  at  the  0.01  (or  even  a  much  smaller)  significance  level.  This  event  and  the 
960729  Lop  Nor  event  illustrate  the  necessity  and  utility  of  regional  event  characterization 
analyses,  particularly  for  cases  in  which  focal  depth  and  Ms:mb  do  not  provide  useful  evidence. 
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Probability  Density 


Figure  32.  Distance-corrected  Pn/Lg  and  Pn/Sn  values,  computed  for  station  NIL  in  the  2-4,  4-6  and  6-8  Hz 
bands,  for  presumed  earthquakes,  two  Lop  Nor  underground  nuclear  explosions  and  the  961110 
Southern  India  event.  The  values  for  the  961110  Southern  India  event  appear  to  be  consistent  with 
other  regional  seismic  activity,  and  considerably  different  than  the  Lop  Nor  nuclear  explosions. 


Figure  33.  Regional  population  analysis  results  for  station  NIL  for  the  961110  Southern  India  event  and  two 
Lop  Nor  nuclear  explosions.  The  event  characterization  parameters  used  are  listed  in  the  legend. 
The  961110  event  is  not  an  outlier  of  the  regional  event  population  at  the  0.01  significance  level. 
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Section  5 

Analysis  of  the  16  August  1997  Event  in  the  Kara  Sea 


5.1  Introduction. 

Based  on  automatic  and  analyst-reviewed  processing  at  the  PIDC,  an  mb  3.9  seismic  event 
occurred  on  16  August  1997  (970816)  at  approximately  02:10:59.9  UTC,  72.65°  latitude  and 
57.35°  longitude  in  the  Kara  Sea  near  Novaya  Zemlya  (NZ).  Figure  34  shows  the  Reviewed  Event 
Bulletin  (REB)  for  this  event,  including  the  associated  arrivals.  The  event  was  approximately 
10.62°  from  NRI,  11.44°  from  SPITS,  16.29°  from  ARU,  16.30°  from  FINES,  20.85°  from  HFS, 
and  21.00°  from  NORES.  ARCES  was  not  operating  at  the  time. 

After  extensive  review  and  use  of  the  Joint  Hypocenter  Determination  (JHD)  program,  using  data 
from  additional  seismic  stations  and  known  event  locations  at  NZ,  the  refined  location  estimate 
remains  offshore  in  the  Kara  Sea,  more  than  100  km  from  the  NZ  test  site,  and  is  consistent  with 
the  original  estimate  in  the  REB  (Israelsson  et  al.,  1997;  Bowers  et  al.,  1997).  It  is  also  near  the 
location  of  a  presumed  earthquake  (mb  4.6)  on  1  August  1986  (Marshall  et  al.,  1989),  as  shown  in 
Figure  35.  Subsequent  analysis,  including  non-IMS  seismic  data,  revised  the  magnitude  estimate 
to  mb  3.3.  A  smaller  (mb  2.5)  event  at  the  same  approximate  location  and  with  similar  seismic 
characteristics  occurred  about  four  hours  after  the  main  970816  event,  indicative  of  a  possible 
aftershock  (Ringdal  et  al.,  1997;  Richards  and  Kim,  1997). 

The  focal  depth  of  the  event  was  constrained  to  the  surface  and  no  long-period  surface  waves  were 
detected.  Hence,  screening  criteria  based  on  depth  and  mb-Ms  could  not  be  applied  and  the  970816 
event  was  not  screened  out,  based  on  insufficient  data.  This  event  was  not  an  outlier  of  the  regional 
earthquake  population  recorded  by  NRI.  The  majority  of  the  regional  events  available  for  such 
comparisons  at  FINES  and  SPITS  occurred  along  the  Mid  Atlantic  Ridge.  Thus,  valid  application 
of  the  outlier  test  could  not  be  made  for  these  stations  due  to  significant  path  differences  for  the 
available  reference  events. 

In  Section  5.2,  all  available  waveforms  for  this  event  are  presented  and  discussed,  including  those 
from  IMS  (NRI,  SPITS,  ARU,  FINES,  HFS  and  NORES)  and  non-IMS  (NORSAR,  KEV,  and  ten 
additional  stations  in  Finland)  stations.  Available  waveforms  from  some  of  these  stations  for 
known  nuclear  explosions  at  NZ,  three  events  near  NZ,  and  regional  earthquakes  are  compared  to 
the  970816  event.  In  Section  5.3,  a  statistical  outlier  test,  using  distance-corrected  high-frequency 
Pn/Sn,  is  used  to  quantify  the  comparisons  to  previous  earthquakes  and  explosions  at  relevant 
regional  stations,  KEV  and  NRI.  In  Section  5.4,  conclusions  regarding  this  event  are  provided. 
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Figure  34.  Reviewed  Event  Bulletin  (REB)  for  the  970816  Kara  Sea  event.  Included  is  a  map  showing  the 
locations  of  the  stations  with  associated  phases  and  the  estimated  location  of  the  event.  Origin  and 
associated  arrival  information  are  listed  below  the  map. 
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5.2  Waveform  Data  for  the  970816  Event. 


Table  9  summarizes  the  available  seismic  data  for  the  970816  event  and  for  other  relevant  nuclear 
explosions  and  earthquakes  used  for  comparison.  Figure  35  depicts  the  locations  of  the  events  near 
or  on  Novaya  Zemlya.  Columns  in  Table  9  include:  (1)  the  stations  (those  with  asterisks  are  IMS 
stations);  (2)  the  distances  of  the  970816  event  to  the  stations;  (3)  the  available  data  for  the  970816 
event  (including  waveforms  and  Pn/Sn  measurements  at  some  stations);  (4)  available  data  for 
previous  NZ  nuclear  explosions;  (5)  available  data  for  other  regional  events  (including  the  860801, 
950613  and  9601 13  events  near  NZ);  (6)  comments  about  the  data  quality  at  the  various  stations; 
and  (7)  a  summary  of  the  regional  population  analysis  results  at  relevant  stations.  The  stations  are 
ordered  by  distance  from  the  970816  event.  The  entries  in  Table  9  that  are  shaded  correspond  to 
cases  for  which  no  applicable  or  valid  comparison  can  be  made. 

Note  that  at  ARU  the  data  for  the  970816  event  had  very  poor  signal-to-noise  ratio  (SNR)  and  Pn 
was  not  detected,  even  after  review  by  an  expert  analyst.  Thus,  no  comparisons  were  made  at  ARU. 

At  KBS  and  NUR  the  data  for  the  NZ  nuclear  explosions  are  all  long-period  (LP).  No  LP  signals 
were  detected  for  the  970816  event.  Thus,  no  comparisons  could  be  made  at  these  stations. 

At  NRI  the  data  for  the  previous  NZ  nuclear  explosions  were  hand-digitized  and  are  not  valid 
above  4  Hz.  Nevertheless,  comparisons  have  been  made  to  the  970816  event  for  the  sake  of 
completeness,  but  with  the  caveat  that  these  comparisons  should  not  be  regarded  as  valid  above  4 
Hz,  the  frequency  range  relevant  to  regional  event  discrimination. 

At  KEV,  NORES  and  NORSAR  there  are  digital  waveforms  for  previous  NZ  underground  nuclear 
explosions  that  are  available  for  valid  comparisons.  There  are  five  additional  NZ  explosions 
recorded  by  NORSAR  (780810,  801011,  821011,  830925,  841025)  that  are  not  listed  in  Table  9 
because  the  records  are  clipped.  The  Sn  signals  for  the  970816  event  at  NORES  and  NORSAR 
have  poor  SNR,  due  to  its  small  magnitude  and  strong  Sn  attenuation  at  far  regional  to  teleseismic 
distances  in  this  region.  Waveform  comparisons  at  NORES  and  NORSAR  are  made,  but  the 
comparisons  shown  below  are  inconclusive,  due  to  poor  Sn-SNR  for  the  970816  event. 

Based  on  the  available  data,  comparisons  of  waveforms  for  the  970816  event  to  those  for  previous 
nuclear  explosions  and  earthquakes  are  made  for  stations  KEV,  NRI,  SPITS,  FINES,  HFS, 
NORES,  and  NORSAR.  Of  these,  the  regional  outlier  analysis,  as  described  in  Section  5.3,  is 
performed  for  stations  KEV  and  NRI. 
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Table  9.  Summary  of  data  for  comparison  to  the  970816  Kara  Sea  event 


Station 

Distance 

(Degrees) 

970816 

Event 

NZ 

Explosions 

Regional 

Earthquakes 

Comment 

Population 

Analysis 

KEV 

9.42 

Waveforms 

Pn/Sn 

821011.841025 
870802,  880507 
881204,  901024 

0 

Good  SNR  for  both 

Pn  and  Sn 

Inconsistent  with 
NZ  EXs  at  0.01 
significance  level 

NRI* 

10.62 

Waveforms 

Pn/Sn 

771009,  801011 
(hdsz  only) 

21  EQs 

hdsz  records  for  EXs 
invalid  above  4  Hz 

Consistent  with 
regional  EQs 

11.44 

Waveforms 

Pn/Sn 

0 

950613 

960113 

Good  SNR  for  both 

Pn  and  Sn 

N/A 

SDF 

11.68 

Waveforms 

0 

0 

No  records  for 
comparison 

N/A 

KBS 

12.41 

Waveforms 

4  (1*  only) 

0 

No  SP  records  for 
comparison 

N/A 

JOF 

13.70 

Waveforms 

0 

0 

No  records  for 
comparison 

N/A 

KJN 

13.71 

Waveforms 

0 

0 

No  records  for 
comparison 

N/A 

SUF 

15.23 

Waveforms 

0 

0 

No  records  for 
comparison 

N/A 

KAF 

15.68 

Waveforms 

0 

0 

No  records  for 
comparison 

N/A 

VAF 

15.94 

Waveforms 

0 

0 

No  records  for 
comparison 

N/A 

KEF 

16.03 

Waveforms 

0 

0 

No  records  for 
comparison 

N/A 

FINES* 

16.30 

Waveforms 

0 

950613 

1800  km  from  FINES 
Strong  Sn  attenuation 

N/A 

PVF 

17.11 

Waveforms 

0 

0 

No  records  for 
comparison 

N/A 

ARU* 

17.15 

Waveforms 
Poor  SNR 

16  (hdsz  only) 

0 

hdsz  records  for  EXs 
invalid  above  4  Hz 

N/A 

NUR 

17.46 

Waveforms 

8  (1*  only) 

0 

No  SP  records  for 
comparison 

N/A 

PKK 

17.93 

Waveforms 

0 

0 

No  records  for 
comparison 

N/A 

HFS* 

20.85 

Waveforms 

0 

950613 
(Poor  SNR) 

Beyond  20  degrees 

No  regional  phases 
Strong  S  attenuation 

N/A 

21.00 

Waveforms 

841025,  870802 
880507,  881204 
901024 

860801 
950613 
(Poor  SNR) 

Beyond  20  degrees 

No  regional  phases 
Strong  S  attenuation 

N/A 

NORSAR 

21.77 

Waveforms 

761020  (sz) 
771009  (sz) 
840826  (sz) 

860801 

Beyond  20  degrees 

No  regional  phases 
Strong  S  attenuation 

N/A 
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Figures  36-43  show  the  available  waveform  data  for  the  970816  event  at  IMS  seismic  stations 
(NRI,  ARU,  SPITS,  FINES,  HFS  and  NORES),  as  well  as  non-IMS  seismic  stations  (NORSAR 
and  11  stations  in  Finland).  All  of  the  waveforms  shown  in  these  figures  have  been  bandpass 
filtered  in  the  4-8  Hz  band. 

Figure  36  shows  that  the  4-8  Hz  Pn  and  Sn  signals  for  the  970816  event  are  comparable  in 
amplitude  on  all  three  components  at  NRI.  Figure  37  shows  that  the  signals  recorded  by  ARU  have 
very  poor  SNR,  particularly  for  Pn. 

Figure  38  shows  that  the  4-8  Hz  Sn  signals  are  smaller  than  the  corresponding  Pn  signals  on  the  sz 
and  bz  channels  at  SPITS,  while  the  4-8  Hz  bandpass  filtered  Sn  signals  are  significantly  larger 
than  the  corresponding  Pn  signals  on  the  bn  and  be  channels  at  SPITS. 

Figure  39  shows  the  4-8  Hz  bandpass  filtered  waveforms  from  FINES.  The  maximum  Sn  signals, 
while  prominent,  are  smaller  than  the  maximum  Pn  signals  on  all  of  the  vertical  channels. 

Figures  40-42  show  4-8  Hz  bandpass  filtered  waveform  segments  from  HFS,  NORES  and 
NORSAR,  respectively,  all  of  which  are  approximately  21  degrees  from  the  location  of  the  970816 
event.  Note  that  only  teleseismic  P  arrivals  were  associated  to  this  event  at  these  stations.  The  P 
waves  appear  impulsive  and  there  are  no  prominent  S  waves  above  the  noise.  These  characteristics 
are  typically  thought  to  be  consistent  with  signals  from  underground  explosions,  but  it  will  be 
shown  that  regional  earthquakes  at  similar  distances  and  locations  also  exhibit  this  behavior. 

Figure  43  shows  the  4-8  Hz  bandpass  filtered  waveforms  from  eleven  Finnish  stations,  including 
KEV.  At  the  present,  we  have  available  data  for  previous  NZ  explosions  only  from  KEV  and  we  do 
not  have  any  earthquake  data  from  these  stations  for  comparison.  However,  these  records  show  that 
for  many  of  these  stations  the  Sn  signals  for  the  970816  event  are  comparable  to  or  larger  than  the 
Pn  signals  (e.g.,  JOF,  KEF,  KEV,  KJN,  SDF,  SUF).  These  waveforms  further  illustrate  that  there 
are  distance  and  station-dependent  variations  in  the  signals. 

Figures  44-56  show  comparisons  of  the  waveforms  for  the  970816  event  to  NZ  nuclear  explosions 
and  regional  earthquakes  recorded  by  KEV,  NRI,  SPITS,  FINES,  HFS,  NORES,  NORSAR.  All  of 
the  waveforms  in  these  figures  are  filtered  in  the  4-8  Hz  passband. 
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970816  Event  (4-8  Hz  Bandpass  Filtered  Waveforms)  -  NRI 


970816  Event  (4-8  Hz  Bandpass  Filtered  Waveforms)  -  ARU 


02:14:00  :16:00  :18:00  :20:00 


Time  (hr:min:sec) 

Figure  37.  Bandpass-filtered  (4-8  Hz)  seismograms  recorded  by  ARU  for  the  970816  event 
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Figure  39.  Bandpass-filtered  (4-8  Hz)  seismograms  recorded  by  FINES  for  the  970816  event. 


Figure  40.  Bandpass-filtered  (4-8  Hz)  waveform  segments  from  HFS  for  the  970816  event. 


970816  Event  (4-8  Hz)  -  Finnish  Stations 


Time  (hr:min:sec) 


Figure  43.  Bandpass-fiitered  (4-8  Hz)  seismograms  recorded  by  eleven  Finnish  stations  for  the  970816  event. 


5.2.1  Waveform  Comparisons  at  KEV. 

Figure  44  shows  vertical-component  (4-8  Hz)  waveforms  at  KEV  for  six  NZ  explosions  (821011, 
841025,  870802,  880507,  881204,  901024)  and  the  970816  event.  The  maximum  Pn  amplitudes 
are  scaled  the  same  for  comparison.  Sn  for  the  970816  event  is  considerably  larger  than  Sn  for  the 
explosions.  Pn/Sn(4-6  Hz)  for  the  970816  event  is  0.90,  less  than  half  the  smallest  value  for  the 
explosions,  which  are  between  1.95  and  3.39.  Figure  45  shows  horizontal-component  (4-8  Hz) 
waveforms  at  KEV  for  three  NZ  explosions  and  the  970816  event.  As  for  the  vertical  components, 
the  970816  event  has  much  smaller  Pn/Sn  values  than  the  NZ  nuclear  explosions. 

Figure  46  shows  unfiltered  and  4-8  Hz  bandpass  filtered  waveforms  for  the  six  nuclear  explosions 
on  all  available  bz  and  sz  channels  at  KEV.  There  were  no  bz  recordings  for  the  821011,  841025 
and  870802  explosions  and  there  were  no  sz  recordings  available  for  the  970816  event.  To  validate 
the  comparisons  of  bz  and  sz  channels,  note  that  the  880507  and  901024  explosions  were  recorded 
on  both  channels  and  their  Pn/Sn(4-8  Hz)  values  are  equivalent  on  both  for  each  event.  (The  sz 
waveform  for  the  881204  event  contains  spurious  high-frequency  content  in  the  4-8  Hz  passband.) 

A  second  significant  point  to  make  regarding  Figure  46  is  that  although  the  unfiltered  waveforms 
for  the  NZ  explosions  have  larger  Sn  than  Pn  amplitudes,  the  Sn  signals  are  significantly  smaller 
than  the  Pn  signals  in  the  4-8  Hz  band.  Thus,  as  many  other  previous  authors  have  noted,  regional 
Pn/Sn  ratios  typically  provide  valid  discrimination  above  3-4  Hz.  Unfiltered  seismograms  or 
passbands  at  lower  frequencies  can  lead  to  misleading  comparisons. 

5.2.2  Waveform  Comparisons  at  NRI. 

Figure  47  shows  the  4-8  Hz  bandpass  filtered  waveforms  at  NRI  for  two  NZ  explosion  signals 
(771009  and  801011),  a  970314  earthquake  and  the  970816  event.  The  970314  earthquake  had  a 
different  azimuth  to  NRI  than  the  970816  event,  but  they  both  occurred  roughly  10  to  11  degrees 
from  NRI  and  have  similar  magnitudes.  The  970314  and  970816  events  have  similar  Pn/Sn  in  the 
4-8  Hz  band.  The  signals  for  the  two  NZ  explosions  were  hand-digitized;  hence,  the  passband 
filtered  signals  exhibit  spurious  spikes  and  should  not  be  considered  valid  in  this  frequency  range. 

5.2.3  Waveform  Comparisons  at  SPITS. 

Figure  48  shows  a  comparison  of  4-8  Hz  filtered  waveforms  for  the  970816  event  at  SPITS  to  those 
for  the  950613  and  960113  events.  Figure  49  shows  the  SPAO/sz  channels  for  the  three  events  in 
the  4-6  and  6-8  Hz  bands.  All  three  events  have  similar  waveforms,  although  Sn  for  the  970816 
event  is  smaller  than  those  for  the  other  two  events.  Note  that  the  970816  event  was  farther  away 
from  SPITS  than  the  other  events;  hence,  relative  attenuation  of  Pn  and  Sn  must  be  considered. 
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Table  10  lists  the  distances  and  Pn/Sn  values  in  the  4-6  and  6-8  Hz  bands  for  these  three  events 
before  and  after  applying  distance  corrections.  The  slopes  of  the  corrections  for  the  4-6  and  6-8  Hz 
bands  are  1.0  and  0.9,  respectively.  Although  the  970816  event  has  the  largest  Pn/Sn  values  before 
applying  corrections,  its  corrected  values  are  not  inconsistent  with  the  950613  and  960113  events 
near  NZ,  given  the  variation  in  the  Pn/Sn  values.  Even  using  conservatively  low  slopes  of  0.6  and 
0.5,  there  is  insufficient  evidence  to  reject  the  970816  event  as  being  consistent  with  the  950613 
and  960113  events. 


Table  10.  Pn/Sn  values  at  SPITS  before  and  after  applying  distance  corrections. 


Date 

Lat 

Lon 

Distance 

(km) 

4-6  Hz 

Before 

6-8  Hz 

Before 

4-6  Hz 

After 

6-8  Hz 

After 

950613 

75.32 

54.85 

1009.32 

1.27 

1.40 

1.88 

2.00 

960113 

75.31 

52.86 

964.30 

0.85 

1.97 

1.31 

2.93 

970816 

72.65 

57.35 

1264.19 

2.31 

2.52 

2.74 

2.93 

5.2.4  Waveform  Comparisons  at  FINES. 

Figure  50  shows  a  comparison  of  bandpass  filtered  waveforms  for  the  970816  event  at  FINES  to 
those  for  the  950613  event  (FLA0  only).  Both  events  exhibit  similar  waveform  characteristics; 
however,  due  to  poor  SNR  for  both  events,  the  comparison  is  inconclusive. 

5.2.5  Waveform  Comparisons  at  HFS. 

Figure  51  shows  a  comparison  of  bandpass  filtered  waveforms  for  the  970816  event  at  HFS  to 
those  for  the  950613  event  (HFSC2/sz  channel  only).  As  for  FINES,  both  events  exhibit  similar 
characteristics  but,  due  to  the  relatively  poor  SNR  for  both  events,  the  comparison  is  inconclusive. 

5.2.6  Waveform  Comparisons  at  NORES. 

Figure  52  shows  a  comparison  of  filtered  waveforms  (4-8  Hz)  from  NORES  for  the  970816  event, 
the  860801  earthquake,  and  six  NZ  nuclear  explosions  (NRBl/sz  channel).  All  of  the  waveforms 
exhibit  very  impulsive  P  waves,  even  those  for  the  860801  earthquake.  The  explosions  exhibit 
negligible  S  waves.  Although  the  S  waves  for  the  860801  earthquake  are  much  smaller  than  the  P 
waves,  they  are  more  prominent  than  the  S  waves  for  the  explosions.  The  S  waves  for  the  970816 
event  are  obscured  by  noise. 
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To  improve  the  SNR,  beams  were  formed  from  the  NORES  array  data  for  the  same  events 
(Figure  53).  As  for  the  NRBl/sz  channel,  the  beams  indicate  that  all  of  the  events  have  very 
impulsive  P  waves  and  very  small  S  waves,  although  the  860801  earthquake  has  noticeably  larger 
S  waves  than  the  explosions.  The  S  waves  for  the  970816  event  appear  to  be  slightly  greater  than 
the  noise.  This  may  indicate  that  the  970801  event  is  more  consistent  with  the  860801  earthquake 
than  the  NZ  nuclear  explosions,  although  the  evidence,  due  to  relatively  poor  SNR,  is  inconclusive. 

A  very  important  point,  however,  is  that  an  earthquake  at  this  distance  from  NORES  (21  degrees 
or  2331  km)  can  have  very  impulsive  P  waves  and  very  small  S  waves,  due  to  strong  attenuation  of 
S  relative  to  P  at  far  regional  to  teleseismic  distances  in  this  region. 

5.2.7  Waveform  Comparisons  at  NORSAR. 

Figure  54  shows  a  comparison  of  the  4-8  Hz  bandpass  filtered  beam  for  the  970816  event  at 
NORSAR  to  those  for  the  860801  earthquake  and  the  761020,  771009,  840826  NZ  nuclear 
explosions.  As  in  Figure  53,  all  of  these  beams  exhibit  very  impulsive  P  waves  and  very  small  S 
waves,  although  the  S  waves  for  the  860801  earthquake  are  much  more  prominent  than  those  for 
the  NZ  explosions.  The  S  waves  for  the  970816  event  are  obscured  by  noise.  Hence,  as  for 
NORES,  this  comparison  is  inconclusive. 

5.2.8  Summary  of  Waveform  Comparisons. 

Figures  55  and  56  illustrate  a  final  comparison  of  waveforms  for  the  970816  event  and  a  961225 
mb  4.1  event  in  the  Jan  Mayen  Island  region.  The  961225  event  was  offshore,  deeper  than  10  km 
at  the  95%  confidence  level  and  with  mb-Ms  <  1.2  at  the  95%  confidence  level.  Hence,  there  is 
strong  evidence  that  the  961225  event  in  the  Jan  Mayen  Island  region  was  an  earthquake.  Figure  55 
shows  that  Sn  attenuates  rapidly  with  distance  and  may  also  depend  on  azimuth.  Figure  56  further 
illustrates  that  it  is  possible  for  an  earthquake  in  this  region  to  have  very  impulsive  P  waves  and 
very  large  high-frequency  Pn/Sn  at  some  stations,  possibly  due  to  combined  effects  of  strong  Sn 
attenuation  along  certain  paths  and  the  shear-wave  radiation  pattern  of  the  focal  mechanism.  Both 
events  have  at  least  one  station  for  which  Pn/Sn  in  the  4-8  Hz  band  is  close  to  one. 

To  summarize,  there  is  evidence  at  KEV  that  the  970816  event  is  inconsistent  with  six  previous  NZ 
nuclear  explosions.  There  is  no  other  evidence  supported  by  these  seismic  data  at  any  of  the 
stations  that  indicates  that  the  970816  event  is  inconsistent  with  a  natural  event. 
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NZ  EXs  and  970816  Event  (4-8  Hz)  -  KEV 
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Figure  44.  Comparison  of  bandpass-filtered  (4-8  Hz)  vertical-component  seismograms  recorded  by  KEV  for 
the  970816  event  and  six  Novaya  Zemlya  nuclear  explosions. 
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Figure  45.  Comparison  of  bandpass-filtered  (4-8  Hz)  horizontal-component  seismograms  recorded  by  KEV 
for  the  970816  event  and  three  Novaya  Zemlva  nuclear  explosions. 
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NZ  Explosion  Waveforms  (Unfiltered  &  4-8  Hz)  -  KEV 
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Figure  46.  Unfiltered  and  bandpass-filtered  (4-8  Hz)  seismograms  (sz  and  bz  channels)  recorded  by  KEV  for 
six  Novaya  Zemlya  explosions. 
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NZ  EXs,  970314  EQ,  970816  Event  (4-8  Hz)  -  NRI 
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Figure  47.  Comparison  of  bandpass-filtered  (4-8  Hz)  seismograms  recorded  by  NRI  for  the  970816  event,  a 
970314  earthquake,  and  two  Novaya  Zemlya  nuclear  explosions  (hand-digitized). 
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950613,  960113,  970816  Events  (4-8  Hz)  -  SPITS 
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Figure  48.  Comparison  of  bandpass-filtered  (4-8  Hz)  seismograms  recorded  by  SPITS  for  the  950613, 960113 
and  970816  events  near  Novaya  Zemlya.  The  waveform  tags  indicate  the  channel  and  event 
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950613,  960113,  970816  NZ  Events  (4-6  &  6-8  Hz)  -  SPITS 
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Figure  49.  Comparison  of  bandpass-filtered  (4-6  and  6-8  Hz)  seismograms  recorded  by  SPAO  for  the  950613, 
960113  and  970816  events  near  Novaya  Zemlva. 
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970816  and  950613  Events  (4-8  Hz)  -  FIAO 
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Figure  51. 
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Comparison  of  bandpass-filtered  (4-8  Hz)  seismograms  recorded  by  FIAO  for  the  970816  and 
950613  events  near  Novaya  Zemlya. 


970816  and  950613  Events  (4-8  Hz)  -  HFSC2 
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Comparison  of  bandpass-filtered  (4-8  Hz)  seismograms  recorded  by  HFSC2  for  the  970816  and 
950613  events  near  Novaya  Zemlya. 


Comparison  to  NZ  EXs  and  860801  NZ  EQ  (4-8  Hz)  -  NORES 


Figure  52.  Comparison  of  filtered  (4-8  Hz)  seismograms  recorded  by  NRB1  for  the  970816  event,  the  860801 
NZ  earthquake  and  the  841025,  870802, 880507,  881204  and  901024  NZ  nuclear  explosions. 
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NZ  EXs,  860801  EQ,  970816  Event  (4-8  Hz  Beam)  -  NORES 
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Figure  53.  Comparison  of  bandpass-filtered  (4-8  Hz)  beams  recorded  by  NORES  for  the  970816  event,  the 
860801  NZ  earthquake  and  the  841025  and  870802  NZ  nuclear  explosions. 


83 


NZ  EXs,  860801  EQ  and  970816  Event  (4-8  Hz  Beam)  -  NORSAR 
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Figure  55.  Bandpass-filtered  (4-8  Hz)  waveforms  recorded  by  NRI  (10.62°),  SPAO  (11.44°),  FIAO  (16.30°), 
HFSC2  (20.85°),  NRA1  (21.00°)  for  the  970816  event. 
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Figure  56.  Bandpass-filtered  (4-8  Hz)  waveforms  recorded  by  SPAO,  ARA0,  HFSC2,  NRA1,  FIAO  for  an  mb 
4.1  earthquake  in  the  Jan  Mayen  Island  region. 


5.3  Regional  Population  Analysis. 


An  outlier  test,  using  distance-corrected  Pn/Sn,  is  applied  here  to  quantify  comparisons  of  the 
970816  event  to  the  NZ  nuclear  explosions  recorded  by  KEV,  and  to  regional  earthquakes  recorded 
by  NRI.  The  reference  events  used  at  NRI  are  all  above  mb  3.5,  to  eliminate  potential 
contamination  from  mining  blasts,  and  do  not  include  events  near  the  NZ  test  site.  Using  the  leave- 
one-out  procedure,  each  reference  event  is  also  tested  as  an  outlier  of  the  remaining  group.  (Fisk  et 
al.,  1993,  1994,  1995,  1996a,  1996b,  describe  the  outlier  procedure  and  previous  applications.) 

5.3.1  KEV. 

Figure  57  shows  the  results  for  KEV,  comparing  the  970816  event  to  six  previous  NZ  nuclear 
explosions.  On  the  left  is  a  plot  of  log  Pn/Sn  values  in  the  4-6  Hz  band  for  the  six  NZ  explosions 
and  the  970816  event.  (Note  that  two  of  the  NZ  explosions  have  nearly  identical  Pn/Sn  values.) 
The  plot  on  the  right  shows  the  results  of  the  outlier  test  using  log  Pn/Sn(4-6  Hz).  Only  one 
discriminant  was  used  since  there  were  only  six  events  in  the  reference  set.  The  distribution  shown 
corresponds  to  that  of  the  likelihood  ratio  obtained  from  the  Pn/Sn  values  for  the  explosions.  The 
vertical  line  corresponds  to  the  threshold  of  the  test  at  0.01  significance  level.  Figure  57  (right) 
shows  that  the  970816  event  is  an  outlier  of  the  NZ  nuclear  explosions  recorded  by  KEV.  This 
indicates  that  the  970816  event  is  inconsistent  with  these  previous  NZ  nuclear  explosions. 


5.3.2  NRI. 

Figure  58  (left)  shows  plots  of  log  Pn/Sn  in  the  2-4,  4-6  and  6-8  Hz  bands  versus  log  distance  for 
21  regional  events  and  the  970816  event.  Figure  58  (right)  shows  the  results  of  the  outlier  test  using 
log  Pn/Sn  in  the  4-6  and  6-8  Hz  band,  indicating  that  the  970816  event  is  clearly  not  an  outlier  of 
the  regional  earthquakes  recorded  by  NRI.  Figure  59  shows  the  locations  of  the  events  that  were 
recorded  by  NRI.  The  970816  event  and  the  970314  earthquake  included  in  the  waveform 
comparison  of  Figure  47  are  indicated.  Since  the  21  regional  events  recorded  by  NRI  have  a 
significant  range  of  distances  and  azimuths,  it  is  important  to  understand  what  effect  this  may  have 
on  the  distribution  of  the  Pn/Sn  values,  particularly  for  events  along  the  Mid  Atlantic  Ridge  versus 
those  at  other  distances  and  azimuths. 

Figure  60  shows  plots  of  Pn/Sn  in  the  4-6,  6-8  and  8-10  Hz  bands  versus  distance  (left)  and 
azimuth  (right)  for  the  events  recorded  by  NRI.  Markers  shown  in  gray  correspond  to  outliers  that 
were  removed  from  the  reference  set  before  the  970816  event  was  tested.  The  data  exhibit 
relatively  little  variation  in  their  Pn/Sn  values  over  a  relatively  wide  range  of  distances  and 
azimuths.  This  suggests  that  these  events  can  be  treated  as  a  single  earthquake  population  at  NRI. 
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Figure  57.  On  the  left  is  a  plot  of  log  Pn/Sn  values  in  the  4-6  Hz  band  at  KEY  for  six  NZ  nuclear  explosions 
and  the  970816  event.  The  results  of  the  outlier  test  using  log  Pn/Sn(4-6  Hz)  is  shown  on  the  right. 
The  970816  event  is  an  outlier  of  the  NZ  explosion  population  at  the  0.01  significance  level. 


Figure  58.  Plots  of  log  Pn/Sn  values  in  the  2-4,  4-6  and  6-8  Hz  bands  versus  log  distance  at  NRI  for  21 
earthquakes  and  the  970816  event  are  shown  on  the  left.  The  results  of  the  regional  outlier 
analysis  at  NRI  using  log  Pn/Sn  in  the  4-6  and  6-8  Hz  bands  are  shown  on  the  right.  The  970816 
event  is  not  an  outlier  of  the  regional  earthquake  population  recorded  by  NRI. 
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Figure  59.  Map  showing  the  locations  of  the  970816  event  and  the  21  reference  events  used  in  the  regional 
population  analysis  at  NRI.  The  overlay  is  a  tectonic  grid  established  by  Oli  Guudmundsson. 
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Figure  60.  Plots  of  log  Pn/Sn  versus  distance  (left)  and  azimuth  (right)  for  the  events  recorded  by  NRI. 


5.4  Conclusions. 


Table  1 1  summarizes  the  comparisons  that  were  made  of  the  970816  event  to  regional  earthquakes 
and  explosions.  This  investigation  of  the  970816  event  indicates  that  its  high-frequency  Pn/Sn 
values  are  inconsistent  with  six  previous  NZ  nuclear  explosions  recorded  by  KEV  and  that  they  are 
consistent  with  21  regional  earthquakes  recorded  by  NRI.  Comparisons  at  NRI  further  indicate  that 
the  970816  event  has  nearly  identical  high-frequency  Pn/Sn  values  as  a  970314  event  at  a  similar 
distance.  The  comparison  at  NRI  to  hand-digitized  data  for  previous  NZ  nuclear  explosions  is  not 
valid  above  4  Hz,  the  frequency  range  relevant  to  regional  event  discrimination. 


Table  11.  Summary  of  conclusions  for  the  970816  Kara  Sea  event. 


Station  S. 

Distance 

(Degrees) 

Comparison 

Comment 

Conclusion  . 

KEV 

9.42 

Waveforms  and  Pn/Sn 
compared  to  six  NZ  EXs 

Good  SNR  for  both 

Pn  and  Sn 

♦  Inconsistent  with  six  NZ  EXs  at  the  0.01 
significance  level  using  log  Pn/Sn  in  the 

4-6  Hz  band 

NRIa 

10.62 

Waveforms  and  Pn/Sn 
compared  to  regional  EQs 

hdsz  records  for  EXs 
are  invalid  above  4  Hz 

•  Consistent  with  21  regional  EQs  using 
log  Pn/Sn  in  the  4-6  and  6-8  Hz  bands 

SPITS* 

11.44 

Waveforms  and  Pn/Sn 
compared  to  two  NZ  EQs 

Good  SNR  for  both 

Pn  and  Sn 

•  Insufficient  evidence  to  conclude  that  the 
970816  event  is  different  than  the  NZ  EQs 

FINES* 

16.30 

Waveforms  compared  to 
950613  NZ  event 

1800  km  from  FINES 
Strong  Sn  attenuation 

•  Inconclusive  evidence  based  on  poor 

SNR  for  the  970816  event 

HFS* 

20.85 

Waveforms  compared  to 
950613  NZ  event 

Beyond  20  degrees 
Strong  S  attenuation 
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a.  Seismic  stations  of  the  International  Monitoring  System  (IMS). 


Comparisons  at  SPITS  to  the  950613  and  960113  events  near  NZ  indicate  similar  waveforms  in  the 
4-6  and  6-8  Hz  bands.  Although  Pn/Sn  in  these  bands  are  higher  for  the  970816  event  than  for  the 
other  two  events  prior  to  correcting  for  distance,  all  three  events  have  more  comparable  values  after 
applying  distance  corrections.  There  is  insufficient  evidence  to  conclude  that  the  970816  event  is 
significantly  different  than  the  950613  and  9601 13  events. 

Comparisons  at  FINES,  HFS,  NORES  and  NORSAR  to  specific  NZ  explosions  and  other  nearby 
events  are  inconclusive  due  to  poor  SNR,  although  there  may  be  some  evidence  using  beams  at 
NORES  that  the  970816  event  is  more  consistent  with  the  860801  presumed  earthquake  near  NZ 
than  with  the  841025  and  870802  NZ  explosions. 
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Examination  of  4-8  Hz  bandpass  filtered  waveforms  at  eleven  Finnish  stations  showed  that  for 
many  of  these  stations  the  Sn  signals  are  comparable  to  or  larger  than  the  Pn  signals  (e.g.,  JOF, 
KEF,  KEV,  KJN,  SDF,  SUF).  These  waveforms  further  illustrate  that  there  are  distance  and  station- 
dependent  variations  in  the  signals.  At  the  present,  however,  there  are  no  available  explosion  or 
earthquake  data  from  these  stations  for  comparison. 

Based  on  the  available  data,  the  seismic  evidence  supports  the  following  conclusions: 

•  The  970816  event  is  inconsistent  with  previous  known  explosions  at  the  NZ  test  site 
recorded  by  KEV,  the  closest  station; 

•  The  970816  event  is  consistent  with  earthquakes  recorded  by  NRI,  the  next  closest  station; 

•  There  is  insufficient  evidence  to  conclude  that  the  970816  event  is  significantly  different 
than  the  950613  and  9601 13  events  near  NZ  recorded  by  SPITS,  the  third  closest  station; 

•  Evidence  at  the  other  stations,  which  are  at  far  regional  to  teleseismic  distances  (e.g., 
1800-2300  km)  is  inconclusive  due  to  poor  signals  and  strong  Sn  attenuation. 

These  conclusions  are  consistent  with  those  obtained  independently  by  Bennett  et  al.  (1997); 
Bowers  et  al.  (1997);  Richards  and  Kim  (1997);  and  Hartse  (1998). 

Since  this  study  was  originally  conducted,  it  was  discovered  that  a  smaller  (mb  2.5)  event  at 
approximately  the  same  location  and  with  similar  seismic  characteristics  occurred  about  four 
hours  after  the  main  970816  Kara  Sea  event,  indicative  of  a  possible  earthquake  aftershock 
(Ringdal  et  al.,  1997;  Richards  and  Kim,  1997). 

Further  analysis  of  the  waveform  data  indicates  that  this  event  does  not  exhibit  any  of  the  expected 
spectral  behavior  expected  for  an  explosion  in  the  water  column  (Baumgardt,  1998;  Herrin,  1998). 
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Appendix  A 

Preliminary  Regional  Distance  Corrections  for  IMS  Seismic  Stations 


Compiled  in  this  appendix  are  scatterplots  of  the  regional  amplitude  ratios,  Pn/Lg  and  Pn/Sn  in  the 
2-4,  4-6  and  6-8  Hz  frequency  bands,  versus  distance,  for  regional  events  recorded  by  43  IMS 
seismic  stations  (26  Primary,  17  Auxiliary)  between  10  September  1997  and  1  October  1997. 
These  stations  recorded  30  or  more  regional  events  with  SNR  for  Pn  greater  than  or  equal  to  2.0. 
There  were  an  additional  7  IMS  seismic  stations  (1  Primary,  6  Auxiliary)  with  30  or  more  events 
that  are  not  represented  here  because  the  estimated  distance  corrections  exhibit  anomalous 
behavior.  The  convention  for  the  station  names  is  that  used  currently  at  the  PIDC.  The  solid  curves 
in  each  plot  represent  the  best  least-square  fits  of  the  regional  distance  dependence.  The  dashed 
curves  correspond  to  the  worldwide-average  distance  dependence. 

Included  with  each  figure  is  a  table  providing  the  number  of  regional  amplitude  measurements,  n, 
the  estimated  standard  deviations  of  the  data  before,  5,  and  after,  sc,  distance  corrections,  and  the 
best-fit  coefficients,  a ,  (3 ,  and  y . 

The  figures  are  ordered  by  Primary  and  then  Auxiliary  seismic  stations.  Within  the  two  sets  of 
stations,  they  are  ordered  alphabetically. 

Note  that  these  distance  corrections  are  in  the  very  preliminary  stages  of  development  and  have  not 
been  validated.  As  such,  they  are  not  recommended  for  operational  use  at  this  time.  In  addition,  the 
distance  corrections  for  many  of  the  stations  require  additional  data  and  treatment  of  subregional 
effects  in  order  to  adequately  improve  and/or  validate  these  results. 


A-l 


ABKT:  Pn/Lg(2-4  Hz)  ABKT:  Pn/Sn(2-4  Hz) 


ABKT:  Pn/Lg(4-6  Hz)  ABKT:  Pn/Sn(4-6  Hz) 


ABKT:  Pn/Lg(6-8  Hz)  ABKT:  Pn/Sn(6-8  Hz) 


Figure  A-l.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  70  regional 
events  recorded  by  ABKT. 
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Figure  A-2.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  90  regional 
events  recorded  by  ARCES. 
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Figure  A-5.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  32  regional 
events  recorded  by  BOSA. 
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Figure  A-6.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  227 
regional  events  recorded  by  CMAR. 
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Figure  A-7.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  128 
regional  events  recorded  by  CPUP. 
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Figure  A-8.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  70  regiona 
events  recorded  by  DBIC. 
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Figure  A-9.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  129 
regional  events  recorded  by  ESDC. 
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Figure  A- 10.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  144 
regional  events  recorded  by  FINES. 
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Figure  A-12.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  216 
regional  events  recorded  by  ILAR. 
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Figure  A- 17.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  52  regional 
events  recorded  by  MNV. 
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Figure  A-18.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  130 
regional  events  recorded  by  NORES. 
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Figure  A-19.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  93  regional 
events  recorded  by  PDAR. 
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Figure  A-21.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  154 
regional  events  recorded  by  PLCA. 


A-22 


STKA:  Pn/Lg(2-4  Hz) 


STKA:  Pn/Sn(2— 4  Hz) 


STKA:  Pn/Lg(4-6  Hz) 


I  ,  ,  .  ,  I 


00 


STKA:  Pn/Sn(4-6  Hz) 


STKA:  Pn/Lg(6-8  Hz) 


STKA:  Pn/Sn(6-8  Hz) 


0  1000  1 

500  2 

- * - — '  o.oi 

000  2500  0 

500  1000  1500  2 

Pn/Lg(2-4  Hz) 


Pn/Lg(4-6  Hz) 


Pn/Lg(6-8  Hz) 


Pn/Sn(2-4  Hz) 


Pn/Sn(4-6  Hz) 


Pn/Sn(6-8  Hz) 


s 

sc 

a 

& 

y 

0.269 

0.242 

0.509 

2.179 

-2.514E-04 

0.088 

0.070 

1.750 

3.389 

-1.029E-03 

0.338 

0.386 

2.480 

4.554 

-1.465E-03 

0.339 

0.296 

-0.484 

0.493 

3.881E-04 

0.158 

0.109 

1.220 

3.121 

-7.251E-04 

0.360 

0.399 

1.743 

3.997 

-1.015E-03 

Figure  A-22.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  37  regional 
events  recorded  by  STKA. 
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Figure  A-27.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  71  regional 
events  recorded  by  AFI. 
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Figure  A-28.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  215 
regional  events  recorded  by  CTA. 
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Figure  A-29.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  179 
regional  events  recorded  by  DAVOS. 
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Figure  A-32.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  111 
regional  events  recorded  by  HFS. 
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Figure  A-33.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  49  regional 
events  recorded  by  HNR. 
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Figure  A-34.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  118 
regional  events  recorded  by  INK. 
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Figure  A-35.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  46  regional 
events  recorded  by  KIEV. 
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Figure  A-36.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  36  regional 
events  recorded  by  KVAR. 
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Figure  A-37.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  112 
regional  events  recorded  by  NIL. 
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Figure  A-39.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  53  regional 
events  recorded  by  PFO. 
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Figure  A-40.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  96  regional 
events  recorded  by  PMG. 
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Figure  A-41.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  71  regional 
events  recorded  by  SNZO. 


A-42 


n 

s 

a 

P 

Y 

Pn/Lg(2-4  Hz) 

30 

0.702 

0.539 

0.711 

4.030 

-6.028E-04 

Pn/Lg(4-6  Hz) 

26 
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0.363 

2.337 
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-1.612E-03 

Pn/Lg(6-8  Hz) 

21 
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3.885 
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-2.592E-03 

Pn/Sn(2-4  Hz) 

32 

0.520 

0.316 

0.891 

3.259 

-4.310E-04 

Pn/Sn(4-6  Hz) 

28 

0.431 

0.306 

1.958 

4.445 

-1.109E-03 

Pn/Sn(6-8  Hz) 

22 

0.263 

0.150 

2.117 

4.065 

-1.184E-03 

Figure  A-42.  Regional  data  and  best  least-square  fits  (solid  curves)  of  the  distance  dependence  for  33  regional 
events  recorded  by  ULN. 
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